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X^yy ' Abstract. We complete our previous investigations concerning the structure and the stability of "isothermal" 

spheres in general relativity. This concerns objects that are described by a linear equation of state, P = qe, so 
that the pressure is proportional to the energy density. In the Newtonian limit g — » 0, this returns the classical 
isothermal equation of state. We specifically consider a self-gravitating radiation {q = 1/3), the core of neutron 
, stars (q = 1/3), and a gas of baryons interacting through a vector meson field (g = 1). Inspired by recent works, 

^ |. we study how the thermodynamical parameters (entropy, temperature, baryon number, mass-energy, etc) scale 

' ' with the size of the object and find unusual behaviours due to the non-extensivity of the system. We compare 

^ , these scaling laws with the area scaling of the black hole entropy. We also determine the domain of validity of these 

scaling laws by calculating the critical radius (for a given central density) above which relativistic stars described 
' by a linear equation of state become dynamically unstable. For photon stars (self-gravitating radiation), we show 

that the criteria of dynamical and thermodynamical stability coincide. Considering finite spheres, we find that 
the mass and entropy present damped oscillations as a function of the central density. We obtain an upper bound 
for the entropy 5* and the mass-energy Af above which there is no equilibrium state. We give the critical value of 
the central density corresponding to the first mass peak, above which the series of equilibria becomes unstable. 
We also determine the deviation from the Stefan-Boltzmann law due to self-gravity and plot the corresponding 
^vgj , caloric curve. It presents a striking spiraling behaviour like the caloric curve of isothermal spheres in Newtonian 

' gravity. We extend our results to d-dimensional spheres and show that the oscillations of mass-versus-central 

density disappear above a critical dimension dcrit{q)- For Newtonian isothermal stars (g — > 0), we recover the 
critical dimension dcrit = 10. For the stiffest stars (g = 1), we find dcrit = 9 and for a self-gravitating radiation 
(g = 1/d) we find dcrit = 9.96404372... very close to 10. Finally, we give simple analytical solutions of relativistic 
^— ^ I isothermal spheres in two-dimensional gravity. Interestingly, unbounded configurations exist for a unique mass 

Mc = cV(8G). 
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1. Introduction M ~ R, which is the exact scaling of the singular isother- 
mal sphere. As a result, the classical thermodynamic 

Self-gravitating systems have a very strange thermody- limit (CTL) for self-gravitating systems corresponds to 

namics due to the attractive long-range, unshielded na- N +oo in such a way that N/V'^^^ is fixed (de Vega 

ture of the gravitational potential (Padmanabhan 1990, & Sanchez 2002). It is then found that the temperature, 

Chavanis 2006c). In particular, they can display negative the energy, and the entropy scale hke T ~ 1, ~ TV, 

specific heats, inequivalence of statistical ensembles, and and S" ~ TV (Chavanis & Rieutord 2003). When quantum 

phase transitions associated with gravitational collapse, effects are taken into account, the mass scales with the 

Furthermore, since these systems are spatially inhomoge- radius as M ~ R~^, which is the exact mass-radius rela- 

neous and their energy non-additive, the usual thermody- tion for nonrelativistic white dwarf stars. As a result, the 

namic limit {N +oo with N/V fixed) is clearly irrele- quantum thermodynamic limit (QTL) for self-gravitating 

vant and must be reconsidered. As a result, the thermo- fermions corresponds to ^ -|-oo in such a way that NV 

dynamical parameters have unusual scalings with respect is fixed. It is then found that the temperature, the energy, 

to the number of particles or system size (see Sect. 7.1 and the entropy scale as T N'^/^, E ^ A^^/^, and S N 

of Chavanis 2006c). For classical self-gravitating isother- (Hertel & Thirring 1971, Chavanis 2002d). These exotic 

mal spheres, the mass scales like the size of the object scahngs (with respect to conventional thermodynamics) 
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come from the long-range action of gravity, which makes 
the system spatially inhomogeneous and renders the en- 
ergy and the entropy non additive. 

On the other hand, in the context of black hole physics, 
Bekenstein (1973) and Hawking (1975) have shown that 
the entropy of a black hole scales as its area S* ~ i?^ and 
that its temperature (measured from an observer at in- 
finity) scales as T ~ 1/R. These results largely remain 
mysterious and are usually considered to reflect a funda- 
mental description of spacetime at the quantum level. In 
order to explain the area scaling of the entropy, researchers 
have invoked holography (Bousso 2002), quantum gravity, 
string theory, entanglement entropy (Srednicki 1993), and 
brick-wall models ('t Hooft 1985). 

In recent works. Banks et al. (2002) and Pesci (2007) 
considered perfect fluids in general relativity with a linear 
equation of state P — qe and found that for the stiffest 
case compatible with causality (g = 1), the entropy scales 
as the area S ^ Rl^ and the temperature as T ~ 1/R. 
Here, these unusual scalings are a consequence of the non- 
extensivity of the system due to the action of gravity. As 
mentioned above, exotic scalings are also encountered in 
Newtonian gravity. However, that purely classical systems 
can exhibit scaling laws analogous to black holes is intrigu- 
ing and deserves further investigation. The study of these 
systems should evidence what, in the area scaling of the 
black hole entropy, is the reflect of a fundamental the- 
ory of quantum gravity and what mainly stems from the 
non-extensive nature of the system. 

In a preceding paper (Chavanis 2002b); denoted Paper 
I, we studied the structure and the stability of spherically 
symmetric relativistic stars with a linear equation of state 
P = qe. These systems are sometimes called "isothermal" 
(by an abuse of language) because in the Newtonian limit 
q ^ 0, they reduce to classical isothermal spheres de- 
scribed by the Emden equation (Chandrasekhar 1972). 
Just as for their Newtonian counterparts, they must be 
enclosed within a box (of radius R) so as to prevent their 
evaporation and make their total mass flnite. For a given 
volume, we showed the existence of a critical mass-energy 
Mc above which there is no equilibrium state. In that 
case the system is expected to collapse and form a black 
hole. Furthermore, considering the series of equilibria, we 
showed that the mass vs central density relation Af (ep) 
presents a series of damped oscillations. The configura- 
tions of hydrostatic equilibrium become unstable above a 
certain value of the central density eo,c corresponding to 
the maximum mass (first peak) . Secondary extrema in the 
series of equilibria correspond to new modes of instability. 
We also numerically observed that the baryon number vs 
central density N{eo) presents extrema at the same loca- 
tions as the mass M{eo), but we were not able to explain 
this observation. 

A first motivation of the present paper is to clarify 
this result. In Sects. [27X112.31 and in Appendix lD.il apply- 
ing the general argument of Weinberg (1972) to a system 
described by a linear equation of state, we show that the 
Oppenheimer-Volkoff equation of hydrostatic equilibrium 



in general relativity can be obtained by extremizing the 
baryon number N[e] at fixed mass-energy M[e]. From the 
condition SN = jiSM (where fi is a Lagrange multiplier), 
it becomes clear that extrema of M(eo) m the series of 
equilibria correspond to extrema of iV(eo), as we observed 
numerically. We also argue in Sect. I2.1l and Appcndix lD.2l 
that only maxima of baryon number N[e\ at fixed mass 
M[e] are dynamically stable. More precisely, the maxi- 
mization of N[e] at fixed M[e] forms a criterion of formal 
nonlinear dynamical stability for a perfect fluid with re- 
spect to the Einstein equations. In Paper I, we considered 
the linear dynamical stability problem and showed that 
the system becomes unstable after the first mass peak 
in the series of equilibria M(eo). In Appendix ID. 21 of the 
present paper, we show that the first mass peak also cor- 
responds to the point at which the hydrostatic configu- 
rations cease to be maxima of the baryon number N[e\ 
(at fixed mass) and become saddle points. Therefore, the 
conditions of linear and nonlinear dynamical stability co- 
incide. Similar results are obtained for barotropic spheres 
described by the Euler-Poisson system in Newtonian grav- 
ity (see Chavanis 2006a). 

Another motivation of our paper, inspired by the works 
of Banks ct al. (2002) and Pesci (2007), is to investigate 
the scaling behaviour of the thcrmodynamical parameters 
(mass-energy M, baryon number N, entropy S, temper- 
ature T, etc) with the system size R. In particular, for 
R +00, it is found in Sects. [Ml and [23] that M i?, 
iV ~ 5 ~ and T ~ i?"FFr. These scaling laws 

were implicit in our preceding paper, but they deserve to 
be emphasised. We complete previous studies in two re- 
spects. In Sees. 12.61 and 12.71 we use the stability criteria 
obtained in Paper I and Appendix [P] to determine the do- 
main of validity of these scaling laws precisely. We show 
that, for a fixed central density eo, the system becomes 
unstable above a critical radius Rc so that the scaling 
laws only hold approximately close to this maximum ra- 
dius. Pure scaling-law profiles, corresponding to singular 
spheres with infinite central energy or to regular spheres 
with very large radii, are dynamically unstable. Secondly, 
we consider general relativistic systems of astrophysical 
interest described by a linear equation of state for which 
all the parameters entering into the scaling laws (includ- 
ing the multiplicative factor) can be calculated explicitly. 
In Sect. [HI we consider a self- gravitating radiation (pho- 
ton star) corresponding to g = 1/3. This problem was first 
studied by Sorkin et al. (1981) in a seminal paper. For this 
system, we note that the entropy is proportional to the 
particle number (S = XN) and the energy is proportional 
to the mass {E = Mc^). Therefore, the maximization of 
N at fixed M (dynamical stability) is equivalent to the 
maximization of S at fixed E (thcrmodynamical stability) . 
We immediately conclude that the conditions of dynami- 
cal and thcrmodynamical stability coincide. We also note 
that the dynamical variational problem 5N = ijl5M can 
be rewritten in the form of a thcrmodynamical variational 
problem 6S = T~^SE where T = c^/(A^) is to be inter- 
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preted as a temperature (we shall see that T corresponds 
to the temperature at infinity Tq given by the Tolman rela- 
tion). We can therefore apply the results of Paper I to that 
system. For a fixed box radius, wc find that there exists a 
maximum mass = Q.2A65. .. MpR/ Lp (where Mp and 
L p are the Planck mass and the Planck length) and a max- 
imum entropy Sc = 0.6217. .. kB{R/ Lp)^^^ above which 
there is no equilibrium state. Moreover, the series of equi- 
libria M(eo) and ^(eo) present damped oscillations and 
become unstable (dynamically and thermodynamically) 
above a critical central density eo,c = 0.439...c''/(G'i?^) 
corresponding to the first mass peak. New modes of in- 
stability appear at each secondary peaks. These stability 
results, which have been proved analytically in Paper I 
and which are developed in the present paper, complete 
and simplify the early analysis of Sorkin et al. (1981). For 
a fixed central density, the photon star becomes unsta- 
ble above a critical radius Rc = 0.663...(c^/G'eo)^^^. For 
large radii, the entropy scales as S* i?^/^. This differs 
from the black hole scaling, but this is consistent with the 
Bekenstein inequality (Bekenstein 1981). We also deter- 
mine the deviation from the Stefan-Boltzmann law due to 
self-gravity and plot the corresponding caloric curve. In 
Sect. 14.11 we find similar results for a gas of completely 
degenerate ultra-relativistic fermions at T = modelling 
the core of neutron stars. In particular, the baryon num- 
ber scales as iV ~ i?^/^ for R — > -|-oo. Wc also consider, 
in Sect. 14.21 a gas of baryons interacting through a vector 
meson field. This model, introduced by Zel'dovich (1962). 
is described by a linear equation of state with q = I. For 
this system the baryon number scales as the area N ^ R^, 
analogously to the black hole entropy. 

A last motivation of our paper is to extend our re- 
sults to a space of dimension d. Such generalization is 
quite common in black hole physics and quantum gravity, 
where is it advocated that extra-dimensions can appear 
at the micro-scales, an idea stemming from the Kaluza- 
Klein theory. In Newtonian gravity, the influence of the 
dimension of space on the laws of physics has only been 
considered recently (Sire & Chavanis 2002, Chavanis & 
Sire 2004, Chavanis 2004, 2006a, 2006b, 2007a). We have 
shown that the structure of the system is highly depen- 
dent on the dimensionality of space and that the problem 
is very rich because it involves several critical dimensions. 
Therefore, it is natural to complete this type of investiga- 
tions. In continuity with our study of classical isothermal 
spheres (Sire & Chavanis 2002), we show in Sect. [5] that 
the oscillations in the mass-central density proflle disap- 
pear above a critical dimension dcritil) depending on the 
index q. Above this dimension, the configurations arc sta- 
ble for any central density, contrary to the case d < dcrit ■ 
This implies that the pure scaling law profiles correspond- 
ing to the singular solution arc now stable. For Newtonian 
isothermal stars (q — > 0) we recover the critical dimension 
dcrit = 10 (Sire & Chavanis 2002). For the stiffest stars 
(g = 1), we find dcrit = 9 and for a self-gravitating radi- 
ation {q = 1/d), we find dcrit = 9.96404372..., very close 
to 10. The oscillations exist for any q £ [0, 1] when d <9, 



and they cease to exist for any q E [0,1] when d > 10. We 
also note that the dimension d = 2 is critical so that the 
results obtained for d > 2 do not pass to the limit d ^ 2. 
In two-dimensional gravity, we obtain in Sect. [6] analytical 
expressions for the density profile of relativistic isothermal 
spheres for any q. Unexpectedly, they have a finite radius 
for q ^ 1 (and the density profile decreases like a Gaussian 
for 9 = 1) contrary to their Newtonian analogue where the 
density decreases as (see, e.g.. Sire & Chavanis 2002). 
Furthermore, they exist at a unique mass Mc = c^/SG. 
Similarly, classical isothermal spheres in d = 2 exist at a 
unique mass Mc = AksT/Gm for a given temperature, or 
equivalently at a unique temperature fesTc = GMm/A for 
a given mass (see, e.g., Chavanis 2007b). Because of that, 
we may expect that these structures are only marginally 
stable. 

2. Relativistic stars 

In this section, we complete the results obtained in Paper 
I concerning the structure and the stability of relativistic 
stars with a linear equation of state. 

2.1. The equations governing equilibrium 

The condition of hydrostatic equilibrium for a spherically 
symmetric perfect fiuid in general relativity is described 
by the Oppenheimer-Volkoff (1939) equation 



1 - 



2GA/(r) 



dP_ 

dr 



P r GM{r 



47rG 



-Pr 



(1) 



where P{r) is the pressure, e(r) is the energy density, and 

1 r 

Mir) = — / e{r)Aiir^ dr (2) 
c Jo 

is the mass-energy contained within the sphere of radius 
r. If R denotes the radius of the configuration, the total 
mass-energy is M = M{R). We also need the baryon num- 
ber N . It is obtained by multiplying the baryon number 
density n{r) by the proper volume element e^^^'^^/'^Airr'^dr 
and integrating over the whole configuration. This leads 
to an expression of the form 



N 



nir) 



1 



2GM(r) 



-1/2 



47rr dr. 



(3) 



If we restrict ourselves to spherically symmetric con- 
figurations and isentropic perturbations, it can be shown 
that the maximization problem 



Max {N[e,n] \ M[e] = M fixed}. 



(4) 



determines stationary solutions of the Einstein equations 
that are dynamically stable (Weinberg 1972). The criti- 
cal points of the baryon number TV at fixed mass AI for 
isentropic perturbations solve the variational problem 



5N - ijl5M = 0, 



(5) 
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where ^ is a Lagrange multiplier enforcing the conser- 
vation of mass. This variational principle leads to the 
Oppenheimer-Volkoff equation (Weinberg 1972); there- 
fore, it determines steady states of the Einstein equations. 
However, these first-order variations tell nothing about the 
stability of the system. Only maxima of N at fixed mass 
M are dynamically stable, so we must consider the sign of 
the second-order variations of N to settle the stability of 
the system. In Appendix[Dl we consider the case of a linear 
equation of state P — qe (see below) so that the baryon 
number is a functional N[e] of the energy density alone. 
In that case, we can evaluate the second-order variations 
of A'^ and study the stability of the system. 

The maximization problem ^ or (|D.ip is similar to 
the minimization of the energy functional yV[p] at fixed 
mass M[p\ = M for a self-gravitating barotropic gas in 
Newtonian gravity. It is known that this minimization 
problem provides a criterion of formal nonlinear dynam- 
ical stability for the Eulcr-Poisson system (see Chavanis 
2006a). Similarly, we expect that the maximization prob- 
lem (j4|) or (jP.ip provides a criterion of formal nonlinear 
dynamical stability for the Einstein equations. 

2.2. The equation of state 

To close the system of Eqs. ([l])-®, we need to specify 
an equation of state relating the pressure P to the en- 
ergy density e. The first law of thermodynamics can be 
expressed as 



™0 



(6) 



where n is the baryon number density and s the entropy 
density in the rest frame. We assume in the following that 
the term Td{s/n) can be neglected. In that case, the first 
law of thermodynamics reduces to 



de ~ — -^dn. 



(7) 



We now assume a "gamma law" equation of state of the 
form 



P ~ qe with g = 7 — 1. 



(8) 



In that case, Eq. ^ can be integrated at once and wc find 
the polytropic relation 



P = A'n'^, 



(9) 



where AT is a constant. Combining Eqs. ([8]) and ([9]), we find 
that the baryon density is related to the energy density by 



1/7 



J/7 



(10) 



The velocity of sound is given by {dP/deY^^c = g^/^c 
so that the principle of causality requires q < 1. In the 
following, we consider < g < 1 (i.e. 1 < 7 < 2). 

There are two situations where the term Td{s/n) can 
be neglected. The first situation is when T = 0. This is the 



situation that prevails in the core of neutron stars where 
the thermal energy is much smaller than the Fermi energy 
so that the neutrons are completely degenerate. The sec- 
ond situation is when the entropy per baryon s/n = A is 
a constant. This is the case in supermassive stars where 
convection keeps the star stirred up and produces a uni- 
form entropy distribution. This is also the case for a gas of 
self-gravitating photons where the pressure is entirely due 
to radiation (see Sect.[3|). For a self-gravitating radiation, 
the Gibbs-Duhem relation 



-P + Ts + jin, 



(11) 



simplifies itself since the chemical potential for photons 
vanishes (/i = 0). In that case, we obtain 

e = -P + Ts. (12) 

Combining the foregoing relations, we find that 

s(r) / \q T{r)\^ 



n(r) 



A 



l + q K 



P{r) = qe{r) = K 



\q T{t) 



1 



K 



(13) 



(14) 



These relations are valid more generally for any system 
with s/n = X constant and fi ^ 0. This is the case for 
example in the cosmology developed by Banks & Fischlcr 
(2001) based on the holographic principle. Central to their 
discussion is a perfect fluid with equation of state P = e 
and T cx s cx e^/^. 



2.3. The general relativistic Emden equation 

Considering the equation of state ([8]), we introduce the 
dimensionless variables ^, ip and M(^) by the relations 
(Chandrasekhar 1972, Chavanis 2002b) 



eoe 



{4^660(1 + 9) 



4 
c q 



1/2 



and 



M(r) 



47reo 



4 
c q 



47rGeo(l + (?) 



3/2 



M{0- 



(15) 



(16) 



In terms of these variables, the Oppenheimer-Volkoff equa- 
tions II])-© can be reduced to the following dimensionless 
forms 



1 



2q MjQ \ dip _ MjQ 

l + q c j dc e 



q^e~ 



(17) 



(18) 



In the Newtonian limit g — > 0, these equations reduce 
to the Emden equation (Chandrasekhar 1942). Therefore, 
Eqs. p7l) - (fT8)) represent the general relativistic equivalent 
of the Emden equation. It is in this sense that relativis- 
tic stars described by a linear equation of state resemble 
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classical isothermal spheres. Just as in Newtonian gravity, 
there exists a singular solution 



e =7?, where Q = -— — .„ . . 
The singular energy density profile is 



(19) 



(20) 



This solution was first found by Klein (1947) and re- 
discovered by numerous researchers including Misner & 
Zapolsky (1964), Chandrasekhar (1972), etc. Considering 
now the regular solutions of Eqs. pT|l - (fT51) . we can always 
suppose that eo represents the energy density at the centre 
of the configuration. Then, Eqs. ([T7])-([TH]) must be solved 
with the boundary conditions ip{0) = i^'{0) = 0. The cor- 
responding solutions must be computed numerically and 
some density profiles are given in Fig. 8 of Paper I for dif- 
ferent values of q. It has to be noted that the asymptotic 
behaviour for ^ — > -|-oo of the regular solutions behave like 
the singular sphere (|^ . 

Since the energy density decreases as r"^ for r — )■ +oo, 
relativistic stars with a linear equation of state have an 
infinite mass like their Newtonian analogues (see, e.g., 
Chavanis 2002a). This reflects the tendency of the sys- 
tem to evaporate. To circumvent this difficulty, we shall 
enclose the system within a box of radius R. Then, the 
solution of Eqs. (|17 p - p^ must be terminated at the nor- 
malized box radius 



|47rGeo(l + g)l'^' 



R. 



(21) 



It can be noted at that point that a is a measure of the 
central density eo for a given box radius R, or a measure 
of the radius R for a given central density eq. The energy 
density contrast is a monotonic function of a given by 



n 



= ^0 _ p4'{a) 

- e{R) 



(22) 



Finally, it is convenient in the analysis to introduce the 
Milne variables (Chandrasekhar 1942) 



The description of the solutions of the general relativistic 
Emden equation in the Milne plane is given in Paper I. 

2.4. Singular solution 

We first give the values of the thermodynamical parame- 
ters corresponding to the singular solution of the general 
relativistic Emden equation. They present exact scaling 
laws that share some analogies with the scaling laws for 
the entropy and temperature of black holes. We stress from 
the beginning, however, that the singular solutions are un- 
stable. The question of the stability of regular isothermal 
spheres is considered in Sect. 12.61 



The singular energy and pressure profiles are 
P(r) _ qQc^ _2 



47rG(l + g) 



(24) 



Using the relation (|10p . we find that the baryon number 
density and the entropy profile are given by 



n(r) = — 



q^Qc^ 



T g+l , 



AnGKil + q) 
From Eq. the temperature profile is 



T(r) = — — A 
\q 



q^Qc^ 



ATTGK{l + q) 



f <?+! , 



(25) 



(26) 



Finally, using Eqs. (106)-(107) of Paper I, we find that the 
functions v{r) and A(r) determining the metric are given 

by 



where A is a constant and 
2q 



V = 



1 



9 



(27) 



(28) 



Comparing Eqs. (|26p and (P7|) . we check explicitly that 
the Tolman relation (see Appendix [B| is satisfied. 

Suppose now that the singular solution is terminated 
by a box with radius R. According to Eqs. and ([M]) . 
the total mass is 



M 



vQc^ 
2G 



R. 



(29) 



Therefore, we find that the mass M scales linearly with 
the radius R. This is the same scaling as for the sin- 
gular isothermal sphere in Newtonian gravity (see, e.g., 
Chavanis 2002a, b) . This is also the scaling entering in the 
Schwarzschild relation (|A.3p . On the other hand, accord- 
ing to Eqs. ([3]) and ([^5]) . the baryon number and the total 
entropy can be written 



(23) N 



S_ _ 

' A " 

x47r 



3g+ 1 



Q— (1-pQ) 



-1/2 



47rGA'(l 



1/7 



3g+l 

R~. 



(30) 



Finally, using the Tolman relation (|B.7[) and Eq. ([^5]) , the 
temperature measured by an observer at infinity is given 

by 



To 



1 + g 
Xq 



K 



q^Qc^ 



47rGA'(l + (?) 



(31) 



The surface temperature T(R) has the same scaling with 
R as To, differing only in the factor (1 —pQY^'^ according 
to the Tolman relation. 
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2.5. Analogy with the black hole entropy 

Some interesting analogies between relativistic stars de- 
scribed by a linear equation of state and black hole ther- 
modynamics have been discussed previously by Banks et 
al. (2002) and Pesci (2007). Let us briefly develop their 
arguments in connection to the present study. 

When P = e, the velocity of sound {dP/deY^^c is equal 
to the velocity of light so that Eq. ([8]) with q = 1 is the 
stiffest equation of state compatible with the principle of 
causality. In that case, we see from Eqs. ([50)) and pTjl 
that the entropy scales like the area S ^ and that the 
temperature scales like the inverse of the radius T ^ 1/R. 
As noticed by Banks et al. (2002), this is similar to the 
scaling of the Bekenstein-Hawking entropy and tempera- 
ture for black holes (see Appendix [A|). We also note from 
Eq. ([25|l that the entropy profile scales as s{r) ~ r~^. 
As discussed by Pesci (2007), this shows that the total 
entropy is due to the contribution of the whole volume. 
Therefore, for relativistic stars with a linear equation of 
state, the area scaling of the entropy is a volume effect 
resulting from the inhoniogeneity of the system. This dif- 
fers from the result of Oppenheim (2002) who considered 
self-gravitating classical systems reaching the above scal- 
ing laws for entropy and temperature while approaching 
the Schwarzschild radius. In this limit, he showed that 
all entropy lies on the surface and the area scaling of the 
entropy is due to external layers. 

More generally, for a linear equation of state of the 

39+1 

form ([8]), the entropy scales as S* ~ i? 9+1 . For < q < 1, 
the exponent is always smaller than 2 (the area scaling 
law) obtained for q = I. This is in agreement with the 
Bekenstein inequality S < Sbh — fcs7r(i?/ip)^ (where 
Lp is the Planck length). Finally, we note that the en- 



ergy scales like E ~ M ^ 



R and that the tempera- 
This leads to - T"^ 



ture scales like T 
implying negative specific heats C = dE/dT < (see 
Appendix \^ for black holes and Sect. 13.51 for the self- 
gravitating radiation). The mass-energy is positive but 
it decreases with the temperature. In Newtonian gravity, 
isothermal spheres also display negative specific heats but 
for a completely different reason (Lynden-Bell & Wood 
1968). According to the Virial theorem 2K + W = 0, 

we have E = K + W K ^NksT (where K is 

the kinetic energy and W the potential energy) leading 
to negative specific heats C = dE/dT = —^Nkg < 0. In 
that case, the energy becomes more and more negative as 
temperature increases. Therefore, the origin of negative 
specific heats for classical isothermal spheres and black 
holes is different. 



2.6. The mass-energy 

We now consider regular isothermal spheres and discuss 
the stability of the system along the series of equilibria 
and the domain of validity of the scaling laws. It is shown 



in Paper I that the normalized mass-energy is related to 
the parameter a by the relation 

2GM pvo{l-quo) 



X 



Rc^ 



1 



(32) 



where mq = u[a) and vq = v{a) are the values of the Milne 
variables at the box radius R. The foregoing relation can 
be rewritten 
2GM 

-^-Xia). (33) 

It defines the series of equilibria parametrized by a going 
from to +00. For a fixed box radius R, we can write 



.1/2 



with 



|47rG'(l +<?) 



c^q 



1/2 



R. 



(34) 



1/2 . 

Therefore, a ~ aeg is a measure of the central density 
so that Eq. ([55)1 determines the relation between the to- 
tal mass M and the central density eq- The correspond- 
ing curve A/(eo) is plotted in Figs. [T] and [H It presents 
damped oscillations around an asymptote at Al ~ Mg = 
X.'iRc^ / {2G) corresponding to the singular solution (P^ . 
For a — > 0, the density is almost uniform and 

2g 2 



Xia) 



(35) 



3(1 + g) 

For a —>■ +00, the density profile tends to the singular 
sphere ([^0]) and 

x{o^) Xs =PQ- (36) 

The mass is maximum for a certain value of the central 
density corresponding to a = ac (first peak). The maxi- 
mum mass is given by 

^ - Xc- (37) 

The values of Uc, Xc and Xs depend on q. For q = 1/3, we 
have ac = 4.7, Xc = 0.493, p = 1/2, Q = 6/7, Xs = 3/7 
and for g = 1, we have ac = 4.05, Xc = 0.544, p = 1, 
Q = 1/2, Xs = 1/2. There is no equilibrium state with 
M > Mc- For M < Mc, the stable configurations corre- 
spond to a < ac, i.e. to sufficiently small central densities. 
The configurations with a > ac in the series of equilibria 
(i.e. those located after the first mass peak) are dynam- 
ically unstable. New modes of instability appear at each 
mass peak. These results are proved analytically in Paper I 
where a linear dynamical stability analysis of box-confined 
systems with a linear equation of state is performed. This 
study is completed in Appendix ID] where it is shown that 
the first mass peak is also the point at which the steady 
states pass from maxima of A^[e] (at fixed mass) to saddle 
points of N[e] (at fixed mass). Therefore, the system loses 
its stability at the first mass peak in agreement with the 
Poincare turning-point argument (Katz 2003, Chavanis 
2006c) applied to the series of equilibria A/(eo). 
For a fixed central density eg, we can write 

1/2 

a = bR, with b = <i ' S> . (38) 



|47rGeo(l 



c4g 
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Fig. 1. Mass as a function of the central density for a 
fixed box radius. This corresponds to the curve x(a). For 
a system evolving at fixed volume, the system becomes 
unstable at the first mass peak. 
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More generally, the mass-radius relation for a fixed central 
density is given by 



2GbM 



x{bR)hR. 



(42) 



This relation is plotted in Fig. [31 The scaling law M ^ R 
is only valid in the limit R +oo. However, as we have 
seen, the system becomes unstable for R > Rc = cuc/b, i.e. 



for a mass AI > Mc 



'RcXc/^G. Therefore, the solu- 



tions exhibiting a pure scaling law profile (like the singular 
isothermal sphere or the regular isothermal spheres with 
R —^ +oo) are dynamically unstable. However, for stable 
configurations close to the critical radius Rc, we observe 
in Fig. [3] that the linear scaling holds approximately, so 
that the results given in Sect. 12.41 are correct in that sense. 
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Fig. 3. Mass as a function of the box radius for a fixed 
central density. This corresponds to the curve axict)- 



Fig. 2. Mass as a function of the central density for a fixed 
box radius showing the secondary peaks. 

Therefore, a = 6i? is a measure of the system size R so 
that the equation 



M = xia) 



c'R 
2G 



(39) 



determines the relation between the total mass and the 
radius. In the limit a — s- 0, the system is homogeneous 
and, using Eq. ([35]), we obtain the usual scaling for an 
extensive system 



M 



3 c2 



(40) 



If we now consider the limit a — > -l-oo corresponding to the 
singular sphere, using Eq. ([3S)) . we get the non-extensive 
scaling 



M 



Xs 



2G 



{R^ +oo). 



(41) 



2. 7. The baryon number 

Using Eq. ([3]) and introducing the dimensionless variables 
defined previously, we find that the baryon number is 
given by 



TV* 

with 



1 



a 1+"! 



e 1+9 



l-p 



M{0 



-1/2 



ed^, (43) 



(44) 



A(a). 



(45) 



4TrGKR^{l + q)_ 
The foregoing relation can be rewritten 
N 

1/2 

For a fixed box radius, using a = aeq , this equation de- 
termines the relation between the baryon number and the 
central density. The corresponding curve is plotted in Fig. 
[H It presents damped oscillations around an asymptote at 



8 
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N ~ Ns = AgN^ corresponding to the singular solution 
(|30|) . According to Eq. ([5]) the extrema of N{a) occur for 
the same values of a as the extrema of M{a) in the se- 
ries of equilibria. This leads to angular points in the curve 
N{M), as shown in Fig. [S] For a ^ 0, corresponding to 
an almost uniform distribution, wc have 



A(a) 



1 ^ 

— Q!<!+1 . 

3 



and for a — 
have 

A(a) ^ A, 



(46) 

-oo corresponding to the singular sphere, we 



Q~{i-xs) 



-1/2 



(47) 



A,. 



The values of ac, A, 
we have ac — 4.7, A, 



(48) 



3q+l 

The baryon number is maximum for a = ac and its value 
is 

Nc 

and As depend on q. For q = 1/3, 
= 0.925, As (8/21)1/4 ^ 0.7856... 
and for g = 1, we have ac = 4.05, Ac = 0.546, A^ = 1/2. 
There is no equilibrium state with N > Nc- For < 
Nc, the stable configurations correspond to a < a^ i.e. 
to sufficiently small central densities. Configurations with 
a > ac in the series of equilibria are dynamically unstable. 
We can note that, for q = 1, the asymptotic values Xs = 
1/2 and A^ — 1/2 coincide. More precisely, the curves 
x(q;) and A(a) turn out to superimpose for a ^ 1. This 
implies that the curve N{AI) becomes linear in the region 
where a ^ 1, as can be seen in Fig. O This behaviour is, 
however, not general and, for q ^ I, the curves xi^^) ^^'^ 
A(a) arc different (see, e.g.. Figs. 18 and 19 of Paper I for 
9=1/3). 




2GM/RC 

Fig. 5. Baryon number as a function of the mass for a 
fixed box radius. Wc have plotted A(q;) as a function of 
x(q;) so that the curve is parametrized by a. For q = 1, 
we note that, asymptotically, A(q;) ~ xio/-) so that the 
relation M ^ N is linear. We have drawn a line y ~ x for 
comparison. This special behaviour arises only for q = I. 

determines the relation between the baryon number and 
the radius. In the limit a 0, the system is homoge- 
neous and, using Eq. (|46p . we have the usual scaling for 
an extensive system 

N^—noR^ (i?^0). (50) 

If we now consider the limit a -l-oo, corresponding 
to the singular sphere, using Eq. P7|) we find the non- 
extensive scaling 
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Fig. 4. Baryon number as a function of the central density 
for a fixed box radius. We have plotted A(a). 

For a fixed central density, using a ~ bR, the equation 



N ^ A(Q!)47r 



q^c^ 



AT:GK{l + q) 



R~, 



(49) 



N = AAtt 



q^c" 



ATiGK{l + q) 



R— 



(R^ -l-oo). 



(51) 



More generally, for a fixed central density, the baryon 
number is expressed function of the radius by 



1 

47r 



^■KGK{l + q) 



q^c^ 



9+1 3q+l 



{bR)^A{bR). 



(52) 



This relation is plotted in Fig. [S] The scaling law TV ~ 

3q+l 

R 9+1 is only valid in the limit R —f +oo. However, as 
we have seen, the system becomes unstable for R > Rc = 
ac/b, i.e. for a baryon number N > Nc = AcN^{Rc). 
Therefore, the solutions exhibiting a pure scaling law pro- 
file (like the singular isothermal sphere or the regular 
isothermal spheres with R —^ +oo) are dynamically un- 
stable. However, for the stable configurations close to the 

3.1 + 1 

critical radius i?c, we note that the scaling N ^ i? 9+1 



holds approximately, so that the results given in Sect 
are correct in that sense. Finally, in Fig. [71 we have plot- 
ted the curve N{M) for a fixed central density. For small 

mass, we have the scaling N M and for large mass, we 

39+1 

have N ~ M^+^. 
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(1981). Our work confirms and extends the results of these 
authors. 



3.1. Equation of state for ultra-relativistic particles 

Let us consider a perfect gas of non-interacting particles 
that can be relativistic. We call /(p) the numerical den- 
sity of particles with impulse p. The particle number, the 
kinetic energy, and the pressure are 



N = V /dp, 



10 



In(bR) 



Eku-. ^vj fe{p)dp, 



Fig. 6. Baryon number as a function of the box radius for ^ = o / JP~j~dp, 
a fixed central density. Wc have plotted a 9+1 A(q;). 



37 dp 

where the kinetic energy of a particle is 
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e(p) ~ mc^ 




1 



(53) 
(54) 
(55) 

(56) 



In the ultra-relativistic limit, the kinetic energy of a par- 
ticle with impulse p reduces to 

e(p) = pc. (57) 
Therefore, the pressure is related to the kinetic energy by 

1 Ekin 



P = 



3 V 



(58) 



In the local form, this leads to the linear equation of state 



10 



P = qe with q ~ —. 
^ ^3 



(59) 



Fig. 7. Baryon number as a function of the mass for a fixed Note that, in a d-dimensional universe (see Sect. [S]), the 



central density. We have plotted a 9+1 A(a) as a function 
of Oix{oi) so that the curve is parametrized by a. 



As a final remark, we note that the density contrast 
(22) is a monotonic function of a (behaving like the in- 
verse of the density profile reported in Fig. 8 of Paper I) 
that could be used to parametrize the series of equilibria 
M{a) and N{a) in place of a. In particular, the critical 
energy density contrast corresponding to ac is TZc = 30.1 
for q = 1 and TZc = 22.4 for q = 1/3. Configurations with a 
density contrast greater than these values are dynamically 
unstable. 



3. Self-gravitating radiation: photon stars 

In this section, we apply the preceding results to a radia- 
tion that is so intense that self-gravity (in the sense of gen- 
eral relativity) must be taken into account. This is some- 
times referred to as "photon stars" (Schmidt & Homann 
2000). This problem was first considered by Sorkin et al. 



corresponding value of this parameter would be 
1 



3.2. The Stefan-Boltzmann law 



(60) 



We first recall some basic elements of thermodynamics 
that apply to a pure radiation modeled as a gas of pho- 
tons. We first ignore gravitational effects. The distribution 
function of a perfect gas of bosons without interaction is 
given by the Bose-Einstein statistics 



/ 



1 



/j3 g,3£-,. _ 1' 



(61) 



where ^ is the chemical potential. For ultra-relativistic 
particles, we have e = pc. On the other hand, for particles 
without rest mass, like photons, the chemical potential 
/i = 0. Therefore, the distribution function of a gas of 
photons is simply 



/ 



1 



1 



/j3 g/3pc _ I ■ 



(62) 
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The particle number, the kinetic energy, and the pressure 
can be computed from Eqs. ((53)) - ([SS]) . Writing these rela- 
tions locally and using the properties of the Bosc integrals, 
we obtain the standard results 



with 



h^^^^^^^ 90' 



(63) 



(64) 



The relation between the energy and the temperature 
is the Stefan-Boltzmann law. Using the Gibbs-Duhem re- 
lation PT|) with /i = 0, we directly obtain the entropy 



(65) 



For a pure radiation, the pressure is related to the photon 
density by 



P = Kri<, 
with 

4 

7 = o' 



K 



he 



90 [87rC(3)4 



1/3- 



(66) 



(67) 



On the other hand, the entropy is related to the photon 
density by 



with A = 



47r4 



90C(3) 



fcf 



(68) 



If we now consider a gas of self-gravitating photons in 
general relativity (photon stars), the preceding relations 
still hold locally. They have the form considered in Sect. 
12.21 with (7 = 1/3. Therefore, a self-gravitating radiation 
can be studied with the theory developed in Sect. We 
note that the singular profiles of Sect. 12.4] scale like 



Pcxeocr^, scxricxr T cx r ^^"^ . (69) 

These are also the asymptotic behaviours of the regular 
profiles for r +oo. The numerical constants in Eqs. 

can be obtained explicitly by using the values of 
-fC, 7, and A given above. 



1/2 



a = \—r^) R- (72) 

For a given box radius, the mass-central density relation 
M(eo) presents damped oscillations around an asymptote 

Ms = \xsMp^, (73) 

with Xs =3/7 corresponding to the singular solution (see 
Fig. [8]) . There is no equilibrium state with a mass larger 
than 

Mc^^XcMp^, (74) 

where Xc = 0.493 corresponding to etc = 4.7. On the other 
hand, the configurations with central density 

£0 > eo,c = T^— ^ , (75) 
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corresponding to a > ac are dynamically and thermody- 
namically unstable (see Sect. 
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Fig. 8. Mass as a function of the central density for a fixed 
box radius. 



3.3. The mass-energy of self-gravitating photons 



We now apply the general results of Sect. l2.6l to the specific 
context of a self-gravitating radiation. The consideration 
of an explicit physical example allows us to determine the 
numerical value of the multiplicative constants appearing 
in the scaling laws. It is convenient to introduce the Planck 
length and the Planck mass. 



Lp = 



Gh 



1/2 



1/2 



(70) 



In terms of these parameters, the relation ([55)1 can be 
rewritten 

M 



1 i? 



(71) 



For a given central density, the mass-radius relation 
M{R) is represented in Fig. [9l For small radii, we have 
the extensive scaling M ~ and for large radii, we get 
the scaling law 



(76) 



MIR 

Mp = 2^^i;' (^^+-)- 
However, configurations with 

^>^^-"^(l6^j ' 

corresponding to a > aj, = 4.7, are dynamically and ther- 
modynamically unstable. This corresponds to 

1/2 



(78) 
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Fig. 9. Mass as a function of the box radius for a fixed 
central density. 

3.4. The entropy of self-gravitating photons 

The total entropy of the gas of photons is 



s(r) 



Using Eq. 
S = XN, 



2GM{r) 
we find that 



-1/2 



47rr^(ir. 



(79) 



(80) 



where N is the number of photons ([3]). Using the results 
of Sect. 12. 7[ we obtain 



S = Al^{a)kB ( ^ 
Lp 



3/2 



with the numerical constant 
1 / 8^3x1/4 



A 



15 



0.672188... 



(81) 



(82) 



For a given box radius, the entropy vs central density 5(eo) 
presents damped oscillations aroimd an asymptote 



R 



3/2 



(83) 



with As = (8/21)^/"', corresponding to the singular solu- 
tion (see Fig. llOp . There is no equilibrium state with an 
entropy greater than 



R 

Sr. = AArkR ( — 

Lp 



3/2 



(84) 



where Aj, = 0.925 corresponding to ac = 4.7. The configu- 
rations with central density eo > eo,c are dynamically and 
thcrmodynamically unstable. Eliminating cq between Eqs. 
((7T|) and ((8T|) . we obtain the entropy vs mass curve S{M). 
This curve is represented in Fig. [11] Since the peaks of 
entropy coincide with the mass peaks in Figs. [8] and [TOl 



this implies that the entropy vs mass curve presents an- 
gular points at each peakf^. However, only the upper part 
of the curve, corresponding to a < ac, is stable. 
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Fig. 10. Entropy as a function of the central density for 
a fixed box radius. 
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Fig. 11. Entropy as a function of mass for a fixed box 
radius. 



For a given central density, the entropy vs radius re- 
lation S{R) is represented in Fig. [121 For small radii, we 
have the extensive scaling S ^ R^ and, for large radii, we 
get the scaling law 



S = AAskB 



R_ 



3/2 



(i?^ +oo). 



(85) 



^ This is similar to the entropy vs energy curve S{E) in 
Newtonian gravity (see Fig. 4 of Chavanis 2002d). 
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However, configurations with R> Rc are dynamically and 
thermodynamically unstable. This corresponds to 



3/4 



(86) 



Finally, eliminating R between Eqs. ([7T|) and (|5T|) . we ob- 
tain the entropy vs mass curve S{M) represented in Fig. 
[T3l For small mass, we have the extensive scaling S ~ M 
and for large mass, we get 



3/2 



(87) 
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Fig. 12. Entropy as a function of the box radius for a fixed 
central density. 
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Fig. 13. Entropy as a function of mass for a fixed central 
density. 



3.5. Thermodynamical stability and temperature 

We have seen in Sect. 12.11 that the maximization of the 
particles number N at fixed mass-energy M provides a 
criterion of formal nonlinear dynamical stability. Since the 
entropy of the self-gravitating radiation is proportional to 
the particle number (5" = XN) and since the energy is 
proportional to mass {E = Mc^), we conclude that the 
criterion of formal nonlinear dynamical stability is equiv- 
alent to the maximization of the entropy at fixed energy: 

Max {S[e] \ E[e] = E fixed}, (88) 

that is to say, to the criterion of thermodynamical stability 
in the microcanonical ensemble. This proves very simply 
that the criteria of dynamical and thermodynamical sta- 
bihty coincided. Now, using S = XN and E = Mc^, the 
variational principle ^ determining the critical points, 
can be rewritten as 

SS = ^5E, (89) 

with T = c^/(A/i). This can be interpreted as the first 
principle of thermodynamics where T is identified with 
the thermodynamical temperature. Now, using Eqs. (|D.7|) 
and we find that 



T = T{R) 1 - 



2G1\I\ 
Rc^ J 



1/2 



(90) 



where T{R) is the surface temperature. Comparing Eq. 
([50)1 with the Tolman relation (|B.7[) . we observe that T 
is equal to the temperature measured by an observer at 
infinity: 

T = To. (91) 

Using Eqs dSSl), ([21]) and (l33|), we find that the 

temperature T is related to the parameter a by the rela- 
tion 

1/2 

(92) 



where 



(«) = 



vl/2 



■[i-x(«)]'/'- 



(93) 



7e(a)i/4 ' 

For a given box radius, Eq. (|92p determines the rela- 
tion r(eo) between the temperature and the central den- 
sity. This curve (see Fig. [H)) presents damped oscillations 
around the temperature of the singular solution corre- 
sponding to = Q^/^il - Xs)i/2 = (96/343)1/"*. There 
is no ec^uilibrium state above a critical temperature 

1 /L \^/^ 



^ For isothermal spheres in Newtonian gravity (Chavanis 
2006a), we have found that the minimization of the energy 
functional W at fixed mass (formal nonlinear dynamical sta- 
bility) is equivalent to the minimization of the Boltzmann free 
energy Fb at fixed mass (thermodynamical stability in the 
canonical ensemble). 
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Fig. 14. Temperature as a function of the central density Entropy as a function of the temperature for 

for a fixed box radius. This curve represents actually the i'^^™^' have plotted A(a) as a function of 

function dipt). 
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Fig. 15. Temperature as a function of the energy for a 
fixed box radius. We have plotted ^(a) as a function of 
X(a)- 



with ~ 0.897 corresponding to the first peak at a* = 
1.47. Since a* < etc, the dynamical and thermodynam- 
ical instability occurs after the first peak of tempera- 
ture H. Eliminating the central density between Eqs. ((7T|) . 



^ We can wonder whether the turning point of mass-energy 
M(eo) is associated with a loss of microcanonical stability and 
the turning point of temperature r(£o) is associated with a loss 
of canonical stability as for Newtonian isothermal spheres; see 
Antonov (1962), Lynden-Bell & Wood (1968), Katz (1978), 
Padmanabhan (1990) and Chavanis (2002a). Note, however, 
that the mass-energy M[e] is a linear functional of the density 
e, like the mass M[p] in Newtonian gravity and contrary to the 
energy i?[p] in Newtonian gravity. Furthermore, we note that 
a self-gravitating radiation becomes unstable above a critical 
temperature or mass-energy, while a classical self-gravitating 
isothermal gas becomes unstable below a critical temperature 



([8T]) and ([92]) , we obtain the temperature vs energy curve 
T{M) and the entropy vs temperature curve S{T). For 
small values of mass or temperature, self-gravity is negli- 
gible and we recover the Stefan-Boltzmann law. For larger 
values, general relativistic effects must be taken into ac- 
count and the curves form a spiral. Figures [15] and 1161 
show the deviation from the Stefan-Boltzmann law due to 
general relativity. For a system evolving at fixed radius, 
the series of equilibria is stable (dynamically and ther- 
modynamically in the microcanonical ensemble) until the 
first turning point of mass or entropy. The caloric curve of 
Fig. [15] is strikingly similar to the caloric curves T{E) of 
Newtonian isothermal spheres (see, e.g., Chavanis 2002a). 

For a given central density, the temperature vs radius 
plot T{R) is represented in Fig. [IT] For small radii, the 
temperature is an intensive variable T = est, and for large 
radii, we have the scaling law 



'SA 



Mpc^ 



R 



1/2 



(95) 



with 6s ~ (96/343)^/^. However, configurations with R > 
Rc are dynamically and thermodynamically unstable (in 
the microcanonical ensemble). This corresponds to 



T<T, 



1 

3A 



(16^)3/4 



a. 



3/2 



„ll/4 



(96) 



Eliminating the radius between Eqs. ([7T|) and ([5^ . we can 
express the temperature as a function of the mass-energy. 
For small masses we have T ~ est and for large masses we 
obtain a law: 



M 
Mp 



1 



2 f MpC^ 



(97) 



or energy. These are important differences that need to be dis- 
cussed further. 
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However, the configurations are unstable for M > Mc or 
T < Tc. Using E = Mc^, the specific heat of the self- 
gravitating radiation (for sufficiently large masses or suf- 
ficiently low temperatures) is given by 



2\3 



dE 



dT 9A2 ° [kBTf 



< 0. 



(98) 



As for black holes (sec Appendix |A|, the specific heat is 
negative. Finally, eliminating the radius between Eqs. (|8ip 
and (j92p . we can express the entropy as a function of the 
temperature. For small entropies we have T ~ est, and for 
large entropies, we obtain a law: 



Mpc^ 



(99) 



-25 



Fixed e„ 
q=1/3 



-15 



ln(M) 



Stable 




M=1/r 



However, the configurations are unstable for 5* > 5c or 
T < Tc. The curves T(M) and S{T) for a fixed central 
density eg are plotted in Figs. [T8land[T9l 



Fig. 18. Temperature as a function of the energy for a 
fixed central density. We have plotted 6{a) / a^l"^ as a func- 
tion of OLx{o)- 
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o 
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Stable 



q=1/3 
Fixed E„ 



Intensive variable 
T=Cst 






In(bR) 



10 



Fig. 17. Temperature as a function of the box radius for 
a fixed central density. This corresponds to the function 
0(a)/ai/2. 
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S=1/T' 
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Fig. 19. Entropy as a function of the temperature for a 
fixed central density. We have plotted ^[ct)Qt'l'^ as a func- 
tion Qid{a)lo^l'^. 



3.6. Comparison with the black hole entropy 



The scaling laws ([75]) . ([55]) . and ([55)1 that we obtained 
for the mass, the entropy, and the temperature of a self- 
gravitating radiation can be compared with the corre- 
sponding scahng laws (|A.3j) . ()A.6p . and (jA.7P for black 
holes. For the self-gravitating radiation, we have 



Sr 



_ 2 /47r3 
~ 3 V315 



1/4 



3/2 



7 ''-^ (ll 



1/2 



(100) 
(101) 



They behave as M - R, Srad R^^^, and Trad R^^^^ 
For the black holes, we have 



M 



1 ,r R 

-Mp . 

2 Lp 



S 



BH 



nkf 



R 



ksTBH = ^Mpc^^. 

47r R 



(102) 



(103) 



They behave as M ^ R, Sbh ^ R^ , an<^ Tbh ^ R^^ ■ We 
have calculated the exact values of the numerical constants 
appearing in the expressions of the mass, entropy, and 
temperature of a self-gravitating radiation. Comparing 
with black holes, the form of the laws are similar but 
the exponents are different. In particular, we note that 
Srad <x- S^^. For large enough radii R where this scaling 
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law applies we have Srad <C Sbh in agreement with the 
Bekenstein (1973) inequality. 

4. Other explicit examples 

In this section, we give other physical examples where the 
general theory developed in Sect. [2] applies. 

4.1. Neutron stars 

In the simplest models, a neutron star can be viewed as an 
ideal gas of self-gravitating fermions. Their density is so 
large that gravitation must be described in terms of gen- 
eral relativity. The distribution function of the neutrons 
is given by the Fcrmi-Dirac statistics 



/ 



1 



h3 e^^-^' + 1 ' 



(104) 



where /i is the chemical potential and e{p) the energy per 
particle. In the completely degenerate limit (T = 0), the 
neutrons have momenta less than a threshold value po 
(Fermi momentum) and their distribution is / = 
There can only be 2 neutrons in a phase space element 
of size on account of the Pauli exclusion principle. 
Furthermore, in the core of neutron stars, the neutrons 
are ultra-rclativistic so that e(p) = pc. According to Sect. 
13.11 they are described by a linear equation of state of the 
form ([8|) with q = 1/3. On the other hand, using Eqs. 
([55]) and ([55)l . the density and the pressure are given by 
n = 87rpg/(3/i3) and P = 27rcpg/(3/i3). Eliminating po 
between these two expressions, we obtain 



P = Kri*, 
with 

4 1 ^3 



1/3 



he. 



(105) 



(106) 



Therefore, the core of neutron stars can be studied with 
the theory developed in Sect. [2]with q = 1/3. The numer- 
ical values of K and 7 are explicitly given by Eq. p06p . 
The singular profiles are 



3c4 



P = 



1/4 



567rG'r2 ' 



1029G3;i3 



-3/2 



(107) 



(108) 



These are also the asymptotic behaviours of the regular 
profiles for r — > +00. On the other hand, for the singu- 
lar solution, the total mass and the total baryon number 
within a sphere of radius R are 



M 



1 



R 



N = BA, 



R 



3/2 



(109) 



(110) 



with xs = 3/7, = (8/21)1/4 ^j^j 
1/4 



1 fA 
^-3U 



0.35408... 



(Ill) 



The description of an ultra-relativistic and completely de- 
generate gas of neutrons in a box (modelling the core 
of a neutron star) is similar to the description of a self- 
gravitating radiation given in Sect. [3l We just need to 
replace the entropy (proportional to the photon number) 
by the neutron number. 

4.2. Baryons interacting via a vector meson field 

Finally, we consider a model introduced by Zel'dovich 
(1962) that provides the stiffcst equation of state P — e 
compatible with the requirements of relativity theory. It 
consists in a gas of baryons interacting through a vec- 
tor meson field (in addition to the gravitational force). 
In the case of neutron stars, the gravitational contrac- 
tion is balanced by the quantum pressure P = Kn^^^ due 
to the Pauli exclusion principle. In the Zel'dovich model, 
the gravitational contraction is balanced by the pressure 
P = Kn? due to the electrostatic repulsion of the baryons. 
For dense objects (n ^ 1), the form of pressure consid- 
ered by Zel'dovich is expected to prevail over the quantum 
pressure. 

Let us first consider a gas of charged baryons with- 
out gravitational interaction. In the Zel'dovich model, the 
(repulsive) potential of interaction between two charges is 
given by 



(112) 



where ^ = |r — r'|, g is the charge of the baryons and fi the 
quantum of mass of the meson (the mass of the meson is 
TO = fj.h/c). Since T = 0, the energy simply represents the 
rest mass energy of the baryons nnihC^ plus the energy of 
interaction 



1 



nnibC 



+00 



it(^)47r^^(i^ = niTLbC^ 



27m2g2 



The pressure is then given by 
d{e/n) 27^2^2 



P 



d{l/n) 



(113) 



(114) 



It can be seen from Eqs. (|113p and (|114p that, in the limit 
of large n, we have 



(9 = 1)- 



(115) 



This represents the most rigid equation of state com- 
patible with relativity theory since the velocity of sound 
{dP/deY^^c is equal to the velocity of light for this value 
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of q. The value g = 1 is clearly an upper bound 0- On the 
other hand, Eq. (|114p can be rewritten 



P = Kri<, 
with 

7 = 2, 



K 



2f;2 



(116) 



(117) 



Therefore, a gas of baryons interacting through a vector 
meson field can be studied with the theory developed in 
Sect. [5] with q = 1. The numerical values of K and 7 are 
explicitly given by Eq. (|117p . The singular profiles are 



e = P = 



WnG 



1/2 



(118) 



(119) 



These are also the asymptotic behaviours of the regular 
profiles for r +00. On the other hand, for the singular 
solution, the total mass and total baryon number within 
a sphere of radius R can be written as 



M 

N = 



1 



R 



R 



(120) 



(121) 



with Xs = 1/2 and — 1/2. Interestingly, these arc the 
same scalings as for the mass and the entropy of a black 
hole (sec Appendix El . Wc have obtained the exact value 
of the prcfactor. 

A more complete description of this system is provided 
by the analysis of Sect. [D The different curves have been 
plotted for q = 1 so that they are directly relevant to the 
Zel'dovich model provided that the numerical values Eq. 
(|117p are incorporated in the discussion. The relation ([55)1 
can be written 



M I , . R 

m; = 2^(")i;' 

with 



(122) 



(123) 



For a given box radius, the mass-central density relation 
M(eo) presents damped oscillations around the mass M^, 



* Note that the above formulae are valid only for fi > l/R 
where R is the system size. For /x = (Coulomb interaction), 
the potential is not screened and the baryons repel each other 
at large distances so that a pressure must be imposed to retain 
the charges. Now, using the Virial theorem 2K + W = 3PV 
and considering = T = 0, we obtain P = e/3 for the elec- 
tromagnetic field. If we take the contribution of the rest mass 
in e into account, we get P < e/3 so that the upper bound is 
q — 1/3 as advocated by Landau & Lifshitz (1960). The model 
of Zel'dovich (1962) shows that this bound can be exceeded for 
a classical vector field with a mass jj,. 



given by Eq. (|120p . corresponding to the singular solution 
(see Figs. [1] and [2]). There is no equilibrium state with a 
mass larger than 



1 R 
Mc = T^XcMp — , 

Z L/p 



(124) 



where Xc = 0.544 corresponding to ac = 4.05. On the 
other hand, the configurations with central density 



eo > eo,. 



iTTG\R) ' 



(125) 



corresponding to a > ac are dynamically unstable. For 
a given central density, the mass-radius relation M{R) is 
represented in Fig.[3l For small radii, we have the extensive 
scaling M ^ R^, and for large radii, we get the scaling law 
(|120p . However, configurations with 



R > R„ = 



1/2 



(126) 



corresponding to a > are dynamically unstable. This 
corresponds to 



M > Mc = Ixcac ( 

2 VsTrG'^eo 



,8 \ 1/2 



(127) 



On the other hand, the relation (1451) can be written 



N = A(a) 



m^c^ / R 



g^h \L 



(128) 



For a given box radius, the baryon number vs central den- 
sity A^(eo) presents damped oscillations (sec Fig.|4|). There 
is no equilibrium state with N larger than 



Nc 



(129) 



where Ac ~ 0.546, corresponding to ac — 4.05. The con- 
figurations with central density eq > £o,c are dynamically 
unstable. Eliminating eg between Eqs. (|122p and (|128p . 
we obtain the baryon number vs mass curve N{M). Since 
the peaks of baryon number coincide with the mass peaks 
in Figs. [HandlH this implies that the baryon number vs 
mass curve presents angular points at each peak a„ (see 
Fig- El)- For a given central density, the baryon number vs 
radius relation N{R) is represented in Fig. [6l For small 
radii, we have the extensive scaling S ^ R^, and for large 
radii, we get the scaling law (|12ip . However, configura- 
tions with R > Rc are dynamically unstable. Eliminating 
R between Eqs. (|122p and (|128p . we obtain the baryon 
number vs mass curve N{M) represented in Fig. [T] For 
small mass, we have the extensive scaling N ^ M and for 
large mass, we get 



N 



^As m^c^ / M 



(130) 
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5. Relativistic stars in d dimensions 

In this section, we extend the theory of relativistic stars 
with a Unear equation of state to d-dimensional spheres. 
This study completes our investigations (Sire & Chavanis 
2002, Chavanis & Sire 2004, Chavanis 2004, 2006a, 2006b, 
2007a) of the dependence of the laws of physics, regarding 
Newtonian gravity, on the dimension of space. It can also 
have applications in the theory of compact objects where 
extra dimensions can appear on the microscale, an idea 
originating from the Kaluza-Klein theory. 

5.1. The equations governing equilibrium 

Let us consider the (d + l)-dimensional, spherically sym- 
metric metric given by 



(131) 



For a perfect fluid, the Einstein equations of general rela- 
tivity governing the hydrostatic equilibrium of a spherical 
distribution of matter are 



dr 



{re-^)-i3-d)e-^ ^{d-2) 



167rG 

{d - l)c' 



dP 

dr 

e^^ dv 1 , , . , 
r dr 



■r^e, (132) 



(133) 



-A 



' (d-l)c4 



(134) 



Combining Eqs. (|132p - p34p . we obtain the d-dimensioual 
generalization of the Oppenhcimcr-Volkoff equation 

aGM{r) \ dP 



„2j.d 2 j fj^j. 



I , , { , , , aGM{r) 8ttG 1 , , 
^(^ + ^) |(^- 2)^ + ^^^P.|, (135) 



with 



M{r) 



Sd 



dr. 



and 



Sd 



27^ 

r{d/2)- 



16n 



(d-l)5/ 



(136) 



(137) 



These equations are only defined for d > 1. For d = 2, we 
have S2 = 27r and a = 8. For c? = 3, we have S3 — Att and 
a = 2. We also recall that the value of the gravity constant 
G and of its dimension G ~ R'^/ {MT"^) changes with the 
dimension of space d. The dimension d = 2 is critical and 
will be treated separately in Sect. [6l In this section, we 
consider that d> 2. 

Using Eqs. p32)) and p34)) . we find that the metric 
functions A(r) and i'[r) satisfy the relations 



,-A(r) ^ ^ _ <^GM (r) 



(138) 



, 1 J leirPr'' 

dv _ (d-l)(d-2)aA/(r)c^ 

dr r ( r'^-^c'^ 1 A 

d-2 \aGM{r) 



(139) 



In the empty space surrounding the star, P = e = 0. 
Therefore, if M = M{R) denotes the total mass-energy of 
the star, Eqs. (fT^ - P^ become for r > i? 



aGM dv 



dr 



d-2 



aGM 



- 1 



The second equation is readily integrated in 
aGM 



v(r) = In 1 



(140) 



(141) 



Substituting the foregoing expressions for A and v in Eq. 
(|13ip . we obtain the {d+ l)-dimensional generalization of 
the Schwarzschild form of the metric outside a star 



This metric is singular at 



dr^ 



1 aGM 



Rs{d). 



'^dn. (142) 



(143) 



This does not mean that spacetime is singular at that ra- 
dius but only that this particular metric is. Indeed, the 
singularity can be removed by a judicious change of coor- 
dinates (see, e.g, Weinberg 1972). Furthermore, it can be 
shown that, for a spherical system in hydrostatic equilib- 
rium, the radius of the star satisfies 



R > 



iGM 



4(d-l) c2 



(144) 



which is the generalization of the Buchdahl theorem in d 
dimensions (Ponce de Leon & Cruz 2000). Therefore, for 
d > 2, the points outside the star always satisfy r > Rs{d) 
so that the Schwarzschild metric is never singular for these 
stars. 

5.2. The general relativistic Emden equation 

Considering the linear equation of state ([8]), we introduce 
the dimcnsionless variables ^, ip, and M{£,) by the relations 



and 



M(r) 



f jd-iyq 
l8^Geo(l + 9) 



1/2 



Sd^o ( [d- l)c'^q 



inGeoil + q) 



M(0. 



(145) 



(146) 



In terms of these variables, Eqs. (|135p and ()136p can be 
reduced to the following dimcnsionless forms 

2q A/(0 1 dip 



1 + q ^-i-^ 
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dM 



(148) 



The Newtonian limit corresponds to q — *■ (see Appendix 
[C]). Taking the derivative of Eq. ((T48)) . we find that 



, _d-l M" 



(149) 



Substituting Eqs. and ([115)) in Eq. pi7|) . we obtain 

a differential equation for the mass profile 



d-l M" {d-2)^+qM' 



M' 



^M(e) 



(150) 



In addition, the metric functions determined by Eqs. (|138p 
and ()133p can be expressed as 



e-^ = l-^ ^-i-V + est, 



1 + q ^''-2 



1 + 9 



(151) 



where the constant is determined by the matching with 
the outer Schwarzschild solution (|14ip at r = i?. 



5.3. Singular solution 



For d > 2, Eqs. ()147|1 and (|148p admit a singular solution 
of the form 



-l/'s _ Q_ 

~ e 



, where Q = 



2(d-2)(l + g) 
(d-2)(l + g)2+4g' 



(152) 



For this singular solution, the metric is of the form 



e^-1 



Aq 



(d-2)(l + g)2' 



(153) 



Coming back to original variables, the singular energy den- 
sity profile is 



Psjr) ^ jd - l)gQe4 1_ 
q 8TrG{l + q) r^' 



(154) 



Considering the thermodynamical variables of Sect, 
we find the scalingsH 

n{r) = 4^ oc r-2/(9+i) , T{r) cx r-^s/l^+i) ^ (155) 
A 



M cx i?' 



d-2 



Tcxi?-29/(9+i). (156) 



The constants of proportionality can be easily obtained 
from the expressions in Sect. 12.41 For q = 0, correspond- 
ing to Newtonian isothermal stars, we obtain the classical 
scalings S - iV, ~ iV and T - 1 with - i?^-^ 
(Chavanis 2004). For q = 1, corresponding to the stiffest 
stars, we find that M ~ R'^''^ and 5' cx iV cx R'^~^ so 
that the entropy scales like the area in any dimension of 
space (Banks et al. 2002). For g < 1, it scales according 
to a power less than the area. On the other hand, the 
temperature scales like T cx 1/i? in any dimension. For 
q = 1/d, corresponding to a self-gravitating radiation or 
to a neutron star in d dimensions, we have the scaling laws 



P cx e cx 7' 



n cx s cx r ''+1 , T en r "i+i . 



M cx R 



d-2 



d(d-l) 

5 CX iV CX 



T cx R-—. 



(157) 
(158) 



Finally, for the sake of completeness, we give the expres- 
sion of the Stefan-Boltzmann law in d dimensions. From 
the Bose-Einstein distribution with a chemical potential 
fi = 0, we find that the pressure, the density of photons, 
and the entropy density are given by 



P 



1 



{d-iy.Sd, 



(kBTr+'Cid+l), 



(159) 



(160) 



(161) 



The equivalent expressions for neutron stars in d dimen- 
sions are given in Chavanis (2007a). 



^ In particular, the entropy scales with the energy E = Mc? 
like S ^ with b = {dq + d-2) /[{q + l){d-2)]. Interestingly, 
this is the same scaling as the one obtained by Kalyana Rama 
(2007) in a cosmological context once the physical size Lh 
of the horizon in his paper is identified with the box size R 
(K.R., private communication). This author studies the phase 
transition, below a critical temperature, from a universe domi- 
nated by highly excited strings to a FRW universe. He then ar- 
gues that the final spacetime configuration (g, d) that emerges 
should maximize the entropy at fixed energy (with the ad- 
ditional condition d > 3). This leads to g = 1 and d = 3, 
providing a possible explanation of why our universe is three- 
dimensional. According to this approach, we have passed from 
a universe with d = 9 dimensions dominated by strings to a 
three-dimensional universe dominated by radiation, then parti- 
cles. Using completely different arguments, we have also found 
that the dimension d = 3 of our universe is special. For ex- 
ample, compact objects like white dwarf stars would be un- 
stable in a universe with d > 4 dimensions (Chavanis 2007a). 
The modifications of the laws of gravity with the dimension of 
space and the special role played by the dimension d = 3 are 
very intriguing. 
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5.4. Asymptotic behaviour 

Considering now the regular solutions of Eqs. (|147p - (|148p . 
we can always suppose that eo represents the energy den- 
sity at the centre of the configuration. Then, Eqs. (|147p - 
(|148[) must be solved with the boundary conditions 



i;iO) = V'(0) = 0. 



(162) 



The corresponding solutions must be computed numeri- 
cally. However, it is possible to determine the asymptotic 
behaviours explicitly. For ^ ^ 0, 



with 

d-2 



qd 



2d 



(163) 



(164) 



where A.{q) is the discriminant 

A{q) = {d- 4)2g2 + 2d{d - I0)q +{d-2){d- 10). (172) 

This function is itself quadratic and its discriminant is 
S = -128(d - l)(d - 10). We have A{q) +oo for q 
±oo. Therefore, for d > 10, A{q) > for all q. For d < 10, 
A{q) = has two roots. Noting that A (0) = {d-2){d-10), 
we conclude that one root is negative and the other is 
positive. On the other hand, noting that A(l) = A{d — 
9){d — 1), we conclude that the positive root is in the 
range ]1, -t-oo[ for d < 9 and in the range [0, 1] for d > 9. 
Summarising: (i) for d < 9, A{q) < for all q G [0, 1] (ii) 
for d > 10, A(g) > for aU q G [0, 1] (iii) for 9 < d < 10, 
then A{q) < for < g < and A{q) > for < g < 1 
(see Figs. [20l and [2T|) where q^,{d) is the positive root of 



and 



b = 



[d{d + 2)q^ + (2^2 ^4d^8)q + d{d- 2)] (2 - d - qd) 



8d2(d + 2)(l + q) 



(165) 



For ^ +00, the asymptotic behaviour of the solution 
of Eqs. (|147[) - (|148p can be obtained by extending the pro- 
cedure developed by Chandrasckhar (1972) in d = 3. We 
introduce a new variable z defined by 



so that z ^ zq = lnQ for ^ +oo. 
variable, Eqs. (fTi7)) - pi5)) become 



(166) 



In terms of this new 




2q M(0 1 dz 
l + q - > - 

d - 2 + (d 4 



J 

2)q M{0 



1 



q 



= 0, 



(167) 



d-3. 



(168) 



dM 
~d^ 

We set z = zo + f with / ^ 1 and linearize the equations. 
We then find that / satisfies the equation 

{d+l)q+_d-l 

^ de 

2{d-2 + 2dq 



Fig. 20. Discriminant A{q) as a function of q for dif- 
ferent values of the dimension of space d (characteristic 
dimensions are described in the text). The mass-central 
density relation M(eo) presents damped oscillations when 
A(q) < and a monotonic behaviour when A{q) > 0. 



q + 1 

{d-2)q^) 



-e 



A{q)-- 
q*{d) 



0, i.e. 

d(10 



d) + v/32(d- l)(10-d) 



[d-Af 



(173) 



(169) 



The further change of variables ^ = e* transforms Eq. 
(|169p to a linear equation with constant coefficients. We 
find 



For d = 9, = 1 and for d = 10, q^ = 0. When A < 0, 
the asymptotic behaviour of Eqs. (|147p - (|148p for ^ -|-oo 
is 



dt^ 



A 



dq + d-2df 2{d-2 + 2dq + {d- 2)q^ 



d-2+dq 
e 2(9+1) 



COS 



2{q+l 



(174) 



1 



dt 



{i + qY 



-/ = o. 

(170) 



Looking for solutions of the form / oc e*^*, we obtain 
-{d-2 + dq)±./A{q) 



k = 



2(g + i) 



(171) 



and it presents damped oscillations around the singu- 
lar sphere. The curve (|174p intersects the singular so- 
lution (|152p infinitely often at points that asymptoti- 
cally increase geometrically in the ratio 1 : exp[27r(l -f 
q)/y/—A{q)]. As shown in Paper I for d = 3 (see also Sect. 
15. 6p . this property is responsible for the damped oscilla- 
tions of the mass-central density relation ^(eo). When 
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Fig. 21. Phase diagram in the {d, q) plane. The line q^,{d) 
separates the region where the mass-central density re- 
lation M(eo) presents damped oscillations from the re- 
gion where the mass-central density relation M(eo) is 
monotonic. The critical dimension corresponding to a self- 
gravitating radiation is obtained by taking the intersection 
between the line (d) and the curve q = 1 /d. 



A > 0, the asymptotic behaviour of Eqs. (|147p - p48p for 

^ —>■ +00 is 



1 



A 



d — 2 + dq-y/A 



(175) 



and it tends to the singular sphere without oscillating. 
In that case, the mass-central density relation M{eo) is 
monotonic (see Sect. 15. 6| ). For a given value of q, the crit- 
ical dimension above which the oscillations disappear is 
such that A{q) =0 leading to 



tiq) 



(176) 



Let us consider specific systems, (i) For classical 
isothermal spheres {q ^ 0), the critical dimension above 
which the mass-central density relation Af(eo) becomes 
monotonic is such that A(0) = (or (d) = 0) cor- 
responding to dcrit = 10. This returns the result ob- 
tained by Sire & Chavanis (2002) who studied the d- 
dimensional Emden equation. As a consequence, for 2 < 
d < dcrit = 10, the caloric curve (3{E) of Newtonian 
isothermal spheres presents a spiraling behaviour around 
the point corresponding to the singular sphere, while 
for d > dcrit — 10, it tends to the point correspond- 
ing to the singular sphere without spiraling (see Sect. 
15. 7p . (ii) For relativistic stars with the stiffest equation 
of state {q — 1), the critical dimension above which 
the mass-central density relation Af (eq) becomes mono- 
tonic is such that A(l) = (or q^{d) = 1) correspond- 
ing to dcrit = 9. Interestingly, this coincides with the 
critical dimension arising in superstring theory that may 
have some connection to the limit case q = 1 (see, e.g., 



Kalyana Rama, 2006). (iii) Finally, for neutron stars or 
for a self-gravitating radiation {q = 1/d), the critical di- 
mension above which the mass-central density relation 
M(eo) becomes monotonic is such that A{l/d) = (or 
(d) = l/d). It is solution of the fourth degree equation 



- lOd^ + d^- 



16 = 0, 



(177) 



leading to dcrit = 9.96404372... very close to d = 10. As a 
consequence, for 2 < d < dcrit = 9.96404372... the mass- 
radius relation of neutron stars M{R) should present a 
spiraling behaviour around the point corresponding to the 
singular sphere (sec Fig. 2 of Meltzcr & Thorne (1966) and 
Fig. 15 of Paper I), while for d > dc„t = 9.96404372..., 
it should tend to the point corresponding to the singular 
sphere without spiraling. A similar property holds for the 
caloric curve of the self-gravitating radiation (in Sect. 15. 7|) . 

5.5. The Milne variables 

As in the Newtonian theory of isothermal spheres 
(Chandrasekhar 1942), it is convenient to introduce the 
Milne variables 

u=^, v = i^'. (178) 



In terms of these variables, the system of Eqs. (|147p - ()148p 
can be reduced to a single first-order differential equation 
(see Paper I for d = 3). For ^ ^ 0, one has 



2a 

and for ^ 



2aC 



Q 
2 ' 



vs = 2. 



(179) 



(180) 



The solution curves in the (u, v) plane in d = 3 for differ- 
ent values of q are given in Fig. 9 of Paper I. According to 
the discussion in Sect. [5^ they tend to the point {us,Vs) 
corresponding to the singular sphere by forming a spiral. 
By contrast, for d > dcrit{q) (defined in Sect. 15.4]) . they 
reach the singular solution {ua,Vs) directly, without spi- 
rahng (see Fig. 



5.6. The mass-density profile 

If the system is enclosed within a box, the solution of Eqs. 
(|147p - (|148p must be terminated at a radius a given by 

'-^P^X'\. (181) 
(d-l)cig J ^ ' 

According to Eq. (|146p the relation between the total mass 
AI of the configuration and the parameter a is 

2q M{a)R'^-^c' 



M 



(182) 



l + q a'^-"^ aG 
Solving for M in Eq. p47p and taking ^ = a, wc get 

qa^e 



M{a) ai)'{a) - nnP- p-^^<^') 



d — 2 -f paijj'{a) 



(183) 
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Fig. 22. The {u,v) plane for isothermal spheres in gen- 
eral relativity. We have considered the case q = 1/d oi a, 
self-gravitating radiation. For 2 < d < dcrit = 9.96404372, 
the curve forms a spiral around the point (us,Vs) corre- 
sponding to the singular sphere. For d > dcrit it reaches 
that point directly, without spiraling. The case ci = 2 is 
studied in Sect. [6l 

where p is defined by Eq. (|28p . Writing uq — u{a) and 
I'o ~ w(q:), wc obtain 



aGM 



PVq{1- quo) 



2- 



pvo 



(184) 



The curve x(a) gives the mass M(eo) as a function of 
the central density for a fixed box radius. It starts from 
X — for a ~ and tends to an asymptotic value Xs = 
pQ/{d - 2) = iq/[{d - 2)(1 + g)^ + Aq], corresponding to 



the singular sphere, as a 
the singular sphere is 



Ms = Xs 



aG ■ 



The mass associated to 



(185) 



The properties of the mass-central density curve x(a) 
can be obtained by a graphical construction in the Milne 
plane. Let us look for the presence of oscillations by find- 
ing possible extrema of xic^)- Taking the derivative with 
respect to a of Eq. (|182[) . using Eqs. (|148|) and (|183p . and 
finally introducing the Milne variables, we get 



dx 
da 



uqvo - {d-2) 



vo{l- quo) 



d- 2- 



pvo 



(186) 



Therefore, the extrema of the curve x(Q!)j determined by 
the condition dx/da = 0, satisfy 



pvo = {d-2) 



1 



uo 



(187) 



The intersections between this curve and the solution 
curve (u, v) in the Milne plane determine the values of 
a for which x(a) is extremum. We easily check that the 
curve (|187p passes by the point {us,Vs) corresponding to 



the singular solution. Therefore, when A(g) < 0, there is 
an infinity of intersections, and the curve x{'^) presents an 
infinity of damped oscillations around the singular sphere 
Xs- This is the case extensively described in d = 3 (see 
Paper I). When A > 0, there is only one intersection (cor- 
responding to the singular sphere), and the curve xi^") 
increases monotonically up to the singular sphere Xs- The 
corresponding curves for a self-gravitating radiation are 
plotted in Fig. [23l If we extend the dynamical stability 



1.2 



0.8 




ln(a) 



Fig. 23. Mass vs central density for a fixed box radius. We 
have considered the case q= 1/d of a self-gravitating radi- 
ation. For 2 < d < dcrit ~ 9.96404372, the curve presents 
damped oscillations around the mass Xs corresponding to 
the singular sphere. The series of equilibria becomes un- 
stable after the first mass peak at (ac,Xc)- For d > dcrit, 
the curve is monotonic. In that case, all the configurations 
of the series of equilibria (with arbitrary central density) 
are stable up to the singular solution with the maximum 
mass x = Xs- The case d = 2 is studied in Sect. [6l 

analysis (see Paper I and Appendix |DJ in a space of d 
dimensions, we find that the point of marginal stability in 
the series of equilibria precisely corresponds to the crite- 
rion (|187p . Therefore, dynamical instability sets it at the 
turning point of mass, as expected. When /S.[q) < 0, the 
series of equilibria is stable until the first mass peak at {ac, 
Mc), and when A{q) > 0, the series of equilibria is stable 
for all values of the central density (including the singu- 
lar sphere that is marginally stable). This implies that, 
for d > dcritil), the pure scaling laws p56p correspond to 
stable configurations contrary to the case d = 3. 

5.7. Self-gravitating radiation 

In this section, we briefly describe the behaviour of the 
curves plotted in Sect. [3] as a function of the dimension 
of space. For sake of generality we give the formulae for 
any q g [0,1], but in the figures, we focus on the self- 
gravitating radiation q = 1/d. 
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Let us first consider a fixed box radius. In that case, 

1/2 

the parameter a oa eg is a measure of the central den- 
sity. The mass-central density relation M(eo) is plotted in 
Fig. [23l For 2 < d < dcrit{q), it presents damped oscil- 
lations around the mass of the singular sphere, and for 
d > dcrit the convergence to the mass of the singular 
sphere is monotonic. The entropy S = XN, where N is 
the number of particles. 



N 



n{r) 



1 



aGM{r) 



1-1/2 



Sdr'^-^dr 



(188) 
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is given as a function of the central density by the relation 



S cx A(a) 



1 



dg+d-2 

a 1+9 "'0 



e 1+9 



1-p 
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-1/2 
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(189) 



For a —t +00 (singular sphere), we have 



A(a) ^ A, = 



dq + d-2 



(190) 



For d < dcrit, the entropy-central density relation ^(eo) 
presents damped oscillations at the same locations as 
M(eo). This implies that the curve S{M) presents some 
peaks at these points. For d > dcrit-, the oscillations in 
S'(eo) and the peaks in S{M) disappear (see Fig.[Ml). The 
temperature 



T = T{R)dl 



aGM 



(191) 



is given as a function of the central density by the relation 

Q,2<?/(<;+i) 



Tex 61(a) = 



7^(a)'^/(9+ 



Ty[l-X(«)] 



1/2 



(192) 



where TZ{a) = e'^'^"' is the density contrast. For a — > -l-oo 
(singular sphere), we have ^(a) ^ 6s = Q~ [1 — Xs)^^^- 
For d < dcrit, the temperature-central density relation 
T(eo) presents damped oscillations at locations different 
from M{eo). This implies that the curve T{M) forms a 
spiral. For d > dcrit, the oscillations in T(eo) and the 
spiral in the caloric curve T{M) disappear (sec Fig. [25|) . 
This is similar to what happens to the caloric curve P{E) 
of a Newtonian isothermal gas for d > dcrit = 10 (see Sire 
& Chavanis 2002). 

Alternatively, if we fix the central density, the param- 
eter a (X i? is a measure of the system size. The curves 
M{R) cx a'^-^xia), S{R) cx a(''«+''-2)/(9+i) A(a), and 
T{R) cx a~'^'^^^'^'^^^6{a) behave like the ones reported in 
Figs. [9l [121 and fTTl (note that, for d > dcritiq), the small 
oscillations arc suppressed). For large R, the scaling laws 
are given by Eqs. ()155|l - (|156p . In particular, the mass be- 
haves like AI ~ XsR'^~^c^ / {aG) for R +00. We again 
emphasise that for d > dcrit the solutions are stable for 
any radius R, contrary to the case d < dcrit where they 
become unstable for R > Rc (see Sect. [2]). 

Finally, we compare the results obtained previously for 
a self-gravitating radiation in general relativity with the 
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26. The solution of the classical Emden equation in 
u, v) plane as a function of the dimension of space. 




Fig. 27. Mass-central density profile for a fixed tempera- 
ture and box radius (or temperature vs density contrast 
for a fixed mass and box radius) as a function of the di- 
mension of space for Newtonian isothermal spheres. 



Fig. 29. The caloric curve as a function of the dimension 
of space for Newtonian isothermal spheres. 
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Fig. 28. Energy vs density contrast for a fixed mass and 
box radius function of the dimension of space for 

Newtonian isothermal spheres. 



results obtained for Newtonian isothermal spheres (Sire 
& Chavanis 2002). In that case, the critical dimension at 
which the oscillations disappear is dcrit — 10. In Fig. 1261 
we plot the solution curve of the Emden equation in the 
{u,v) plane. In Figs. [27l and[28l we plot the inverse tem- 
perature 77 and the energy A as a function of the parameter 
ct = {SdGPmpoY^^R, which is a measure of the density 
contrast. In Fig. [511 we plot the caloric curve /?(£'). For 
d = 2, the caloric curve tends to a plateau with temper- 
ature Tc = GMm/ AkB and energy i5 — *■ — 00. This corre- 
sponds to the formation of a Dirac peak as T ^ Tc. For 
2 < d < 10, the caloric curve forms a spiral. For d > 10, 
the spiral shrinks to a point. These results are strikingly 



similar to those obtained for a self-gravitating radiation 
in general relativity. 

5.8. Particular values of the parameters 

In this section, we regroup the values of the different pa- 
rameters defined in the text for different values of q and 
different dimensions of space. 

For q = 1 (stiffest stars) , we have p = 1, Q ^ 9s = , 
and Xs = As = In particular, for d = 3, we have p = 
1, Q = 9s ~ ^ and Xs = ~ \- Furthermore, = 4.05, 
Xc ^ 0.544, Ac = 0.546, a, = 1.30, and 6*, = 0.676. On 
the other hand, for dcrit = 9, we have p = 1, Q = 0s = |, 
and Xs As = i. 

For q = 1/d (self-gravitating radiation), we have 
P = dTT' Q = 2rf(rf - 2)(d + l)/[(d -2){d+ 1)2 + 4d], 
Xs = pQ/{d - 2), 9s = Qi/(<i+i)(i _ ^^)i/2^ and As - 
_d±i_Qd/{d+i)(^^ _ ^^)-i/2 In particular, for d = 3, we 

have p = I, Q = 1^ Xs ^ I Os = (^)l/^ and 
As = (|-)^/*. Furthermore, = 4.7, Xc = 0.493, 
Ac = 0.925, a, = 1.47 and 9^ = 0.897. On the other 
hand, for P{dcrit) = 10 (closest upper integer value of the 
critical dimension), we have p = Q — Xs — iffg: 

n _ f 110\l/11^1331U/2 and A - ii/' 110 UO/ll/ 1386 Nl/2 



6. Two dimensional gravity 

In this section, we focus on the dimension d = 2, which 
presents peculiar features and where analytical results can 
be obtained. 
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6.1. Self-confined solutions 



For d ~ 2, the Oppcnhcimcr-VolkofT equations (|135p - (jl36p 
reduce to 

SGMirndP__S^^^^^^^^^ (193) 



1 - 
with 

M{r) 



2tt 



C .10 



dr 



erdr. 



(194) 



From Eqs. (fTSS]) and (fT39l) . we find that the functions 
defining the metric are given by 

167rGPr 

(195) 



_^ _ ^ 8GM{r) dv 
(9 ' dr 



1 



8GA/(r) 



In the empty space surrounding the star, P = e = 0. 
Therefore, for r > R, we have e^^ = 1 — 8GM/c^ and 
u = est. The first relation requires M < Mc = c^/SG. 
In fact, we find in the fohowing examples (valid for a lin- 
ear equation of state) that steady solutions exist only for 
M = Mc = c^/8G. We can wonder whether this result is 
general. 

For a linear equation of state, defining 

1/2 



eoe 



4 

c q 



and 



9 



8nGeo{l + q) 



(196) 



(197) 



(198) 



M(r) - 

the generalized Emden equation takes the form 

i-^M(oj5 = .ee-^ 

1 + q j d^ 

dM 

We note that there is no Newtonian limit (g — > 0) in d = 2 
(see Appendix [C|. The differential equation for the mass 
profile is 

1 M" _ qM' 



(199) 



1 



1+9 



(200) 



This equation can be easily integrated once to yield 



M' 



1 - 



2q 



1 



M 



(201) 



where is a positive constant. This can again be inte- 
grated easily to yield 

1-9 



1 



2q 



1 



■M 



q 



2c2 



(202) 



where B and K (related to Aq) are some constants. Here, 
we assume q ^ 1 (the case q = 1 will be treated sepa- 
rately). Because M = at ^ = 0, we find that B = \. 
Therefore, the mass profile is given by 

1 



M(0 = 



2q 



1 - (1 - K'^i^) 



(203) 



Using Eq. (|199p , we find that the normalized density pro- 
file is given by 



1 



(1 - K'^fy 



(204) 



Thus, we find that the density vanishes at a finite distance 
^0 = ^/K. Using furthermore the fact that ^/'(O) = 0, we 
obtain 

2(1 + 9)1'/' 

q{^-q). 

Therefore, the mass and density profiles can be written!^ 



(205) 



M(0 



2q 



1 - (e/Co)' 



[1 - (e/Co)' 



Returning to original variables, we find that 



M(0 = - 



1 - 1 - (r/ro)' 



E(r) = 6o [l-(r/ro)2]-' 

„4 1 1/2 



47rGeo(l-g). 
The metric is explicitly given by 

= [1 - (r/ro)2] , 

2<7 



1-9 



In [1 - [r/roY 



(206) 
(207) 

(208) 
(209) 
(210) 

(211) 
(212) 



This defines a family of solutions parametrized by the cen- 
tral density. These solutions can have different radii y'o, 
but they all have the same mass M — M{ro) = Mc given 

by 



Me 



c 

8G' 



(213) 



For rp 0, the density profile tends to a Dirac peak. In 
two dimensions, the case of photon stars (self-gravitating 
radiation) and neutron stars corresponds to g = 1/2. On 
the other hand, for the stiffest equation of state corre- 
sponding to g = 1 , Eq. (|20ip simplifies in 



M' 

— = Ao{l - M). 

This can be solved easily to yield 

M{0^l-e'^'-^\ 



(214) 



(215) 



® It is amusing to note that the form of the density profile is 
similar to a "Tsallis distribution" with indexp = 2g/(l+g). For 
g = p = 1, we obtain a "Boltzmann distribution" (|217|l . These 
analogies with generalized thermodynamics (Tsallis 1988) are, 
of course, purely coincidental. They show that the "Tsallis dis- 
tribution" can arise in very different contexts that are not nec- 
essarily related to thermodynamics. 
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where we have used M(0) = 0. The density profile is the 
Gaussian 



Using -0(0) = 0, we get K"^ 1/2 so that 

Returning to original variables, we find that 
e(r) = eoe-('-/^)', 



M{r) 



where 



L = 



8G 



1 



1/2 



(216) 

(217) 
(218) 

(219) 
(220) 

(221) 



STrGeo^ 

is a typical lengthscale. The metric is explicitly given by 
e-"" = e-^'''^^\ (222) 



Again, we find a family of solutions parametrized by the 
central density. It is noteworthy that the total mass M = 
M(+(X)) = Mc is the same for all these configurations and 
is again given by Eq. ()213p . For eq +oo, the density 
profile tends to a Dirac peak. 

These results have to be contrasted from their 
Newtonian counterpart where the density decreases as 
for ^ +00 (see, e.g.. Sire & Chavanis 2002). However, 



the universal mass Mr 



V(8G) seems to be the gen- 



eral relativistic equivalent of the critical temperature T^. = 
GMm/{AkB) or the critical mass = AkBT/{Gm) in 
2D Newtonian gravity for isothermal spheres. Indeed, for 
a fixed temperature, unbounded two-dimensional, self- 
gravitating isothermal spheres have a unique mass Mc 
(see, e.g., Chavanis 2007b). 

Using Eq. (|10p . the baryon number is given by 



TV 



+ 00 



(223) 



with n{r) = (g/i^)i/(9+i)e(r)i/(9+i) . From the analytical 
expressions Eqs. pIS)) (PTT|) or ([^1^1) (g^H), we note that 
^^^~)gA(r)/2 _ ^j^g siai. For I, using Eq. (|210p we 
explicitly obtain 



N= 4, 



(224) 



4G(1 ~q)'' ■ 

The baryon number diverges for eo — > 0, i.e. rg +oo, 
suggesting that the system tends to evaporate (recall that 
stable states tend to maximize N at fixed mass M) Q. For 
q — \, the baryon number diverges, whatever the central 
density. 



^ This result is different from the case of 2D Newtonian grav- 
ity where the Boltzmann free energy Fb[p\ happens to be in- 
dependent of the central density po parametrizing the family 
of isothermal solutions at M = Mc or T = Tc (see Chavanis 
2007b). 



6.2. Box-confined solutions 

If the system is confined within a box of radius R, the 
previous results remain valid for r < R. If vq > R, the 
system is confined by the wall {p{R) ^ 0), and if tq < R, 
the system is self-confined (p(ro) = 0). Let us consider 
here box-confined configurations (rg > R). For q ^ 1, the 
mass-central energy relation for fixed R is given by 



^^-8G 



l-ll-iR/roY 



1/2 



47rGeo(l - q) 



(225) 



(226) 



The mass-central density (for a fixed box radius) is plotted 
in Fig. [3D] and the density profile is plotted in Fig. EH For 
q ~ 1, we have 



M 



-i 

8G V 



1 



-{R/LY 



where 



L = 



1/2 



(227) 



(228) 



The mass-central density (for a fixed box radius) is plotted 
in Fig. [32] and the density profile is plotted in Fig. [33] 
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Fig. 30. For box-confined systems, this figure represents 
the mass as a function of the central density (proportional 
to {R/rof) for q = 1/2. 



6.3. The Milne variables 



The Milne variables are defined by Eqs. (|178p . Using the 
analytical solution (|207p . we find for q^ 1 that 



« = -(l-(C/eo)')^: 
1 



1 - (e/?o)= 



(229) 
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Fig. 31. Density energy profile (in units of /2tiGR^) 
for different values of the central density (proportional to 
{R/roY) for q = 1/2. The profile is self-confined if ro < i? 
corresponding to eo > c^/[47rG(l — q)R^]- For ro ^ or 
Co — > +00, the density profile tends to a Dirac peak. 




R/L 



Fig. 32. For box-confined systems, this figure represents 
the mass as a function of the central density (proportional 
to {R/Lf) for g = 1. 

Eliminating ^ between these two expressions, we obtain 



1 



(230) 



The {u,v) curve is parametrized by ^. For ^ = 0, we have 
(u, v) = (1/g, 0) and for ^ (_o, we have (u, v) = (0, +00). 
For q = 1, using the analytical solution (|217p . we find that 

u = e^«'/^ v^C^. (231) 

Eliminating ^ between these two expressions, we obtain 

w = -21nu. (232) 

The (u, v) curve is parametrized by ^. For ^ = 0, we have 
{u, v) = (1,0) and for ^ — > -l-oo, we have (u, v) = (0, -l-oo). 



Fig. 33. Density energy profile (in units of /SttGR^) 
for different values of the central density (proportional to 
[R/LY) for q = 1. For L ^ or eo ^ +00, the density 
profile tends to a Dirac peak. 



The solution curve in the (m, v) plane is represented in 
Fig. [34l If the system is confined within a box of radius 




Fig. 34. The (w, v) curve of "isothermal" spheres in two- 
dimensional gravity for q= 1/2 and q= 1- 

R, the previous results remain valid for ^ < a < Co where 



87rGeo(l + g)1 



c^q 



R. 



(233) 



From the criterion (|187p , we can extrapolate that the con- 
dition of marginal stability in d = 2 corresponds to uq = 0. 
Therefore, the solutions that arc confined by a box are sta- 
ble, since uq = u{a) > 0, while the self-confined solutions 
are marginally stable since uo = u{£,o) = 0. 
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7. Conclusion 

In this paper, we have carried out a thorough analysis of 
the structure and stability of relativistic stars (and self- 
gravitating radiation) described by a linear equation of 
state. In order to prevent evaporation, we have placed 
these objects in a cavity. We have studied the stability of 
the steady states (i) by linearizing the Einstein equations 
around the stationary solution and studying the sign of 
the squared pulsation (see Paper I) and (ii) by consid- 
ering the sign of the second-order variations of the baryon 
number (see Appendix [D|). Wc have found analytically 
that the instability occurs precisely at the first mass peak 
when we vary the central density for a fixed box radius. 
We have applied these results to explicit examples: self- 
gravitating radiation, the core of neutron stars, stiffest 
stars, Zel'dovich model etc. We have determined their do- 
mains of stability precisely and found some upper boimds 
on the thcrmodynamical parameters like the entropy. Wc 
have obtained scaling laws (with a precise determination 
of the prcfactor) that wc compared with the scaling laws of 
black holes. We have stressed, however, that pure scaling 
laws correspond to unstable configurations and that they 
can hold only approximately close to the critical radius 
i?c- We have generalized our results in d dimensions and 
found two critical dimensions. In d = 2, steady state solu- 
tions exist for a unique value of the mass Mc = (? j (SG) 
and they can be obtained analytically. For this mass, we 
get an infinite family of solutions parametrized by the cen- 
tral density. The density profiles vanish at a finite radius 
for g ^ 1 and decrease as a Gaussian for g = 1. These 
configurations arc probably marginally stable. They may 
evolve dynamically cither toward a completely spread pro- 
file (evaporation) or toward a Dirac peak (black hole) with 
mass Mc (collapse). On the other hand, for d > dcrit{q) 
(where dcritio) is a nontrivial dimension depending on q), 
the mass vs central density profile no longer displays oscil- 
lations, and all the configurations of the series of equilib- 
ria are stable whatever their central density. In that case, 
the pure scaling laws are meaningful. The physical impli- 
cation of this result needs further investigation. For the 
self-gravitating radiation, we have found that the critical 
dimension has the non-integer value dcrit = 9.96404372... 
This is an interesting (and intriguing) result because it 
results solely from the combination of the Einstein equa- 
tions and the Stefan-Boltzmann law, so it has a relatively 
fundamental origin. 

We have also shown that the structure of relativistic 
stars with a linear equation of state in general relativity 
is strikingly similar to the structure of isothermal stars 
in Newtonian gravity. Basically, the analogy stems from 
the rescmblence between the general relativistic Emdcn 
equation and its Newtonian counterpart and the fact that 
they coincide for g — 0. On the other hand, stable rela- 
tivistic stars with a linear equation of state maximize the 
baryon number N[e] at fixed mass-energy M[e\ (formal 
ironlinear dynamical stability for the Einstein equations). 
Similarly, isothermal stars in Newtonian gravity minimize 



the energy functional yV[p] at fixed mass M[p\ (formal 
nonlinear dynamical stability for the Euler-Poisson sys- 
tem) or minimize the Boltzmann free energy functional 
Fb[p\ at fixed mass M[p\ (thcrmodynamical stability in 
the canonical ensemble). As a result, the curves of Figs. 
[T] and [?] are respectively similar to (i) the mass-versus- 
central density M(po) at fixed temperature and volume 
and (ii) the energy functional yV(po) or the Boltzmann 
free energy _F'B(po)-versus-central density at fixed tem- 
perature and volume for Newtonian isothermal spheres 
(Chavanis 2002d). On the other hand, by interpreting the 
mass M[e] as an energy and the baryon number N[e] as 
an entropy, the criterion of formal nonlinear dynamical 
stability in general relativity is equivalent to a criterion 
of thcrmodynamical stability in the microcanonical en- 
semble. Therefore, the curves of Figs. HI [TOl HH and [15] 
are respectively similar to (i) the energy E{TZ) and the 
entropy S'(7?.)-vcrsus-dcnsity contrast at fixed mass and 
volume, (ii) the entropy S'(i?)-versus-energy at fixed mass 
and volume, and (iii) the caloric curve [3{E) at fixed mass 
and volume for Newtonian isothermal spheres (Chavanis 
2002d). Therefore, depending on the interpretation, rela- 
tivistic stars with a linear equation of state share analogies 
with classical isothermal spheres in both microcanonical 
and canonical ensembles. These interesting analogies are 
intriguing and deserve further investigation. 

Note added: Coincidentally, two other authors 
V. Vaganov [arXiv : 0707 . 0864 and J. Hammersley 
arXiv : 0707 . 0961 have, in two recent papers, indepen- 
dently carried out studies of the sclf-gravitating radiation 
in d dimensions (in an asymptotically anti-de Sitter space 
A < 0) related to the one developed in Sect.[5l These au- 
thors find that the oscillations in the mass-central density 
relation disappear above a critical dimension. On the basis 
of numerical calculations, J. Hammersley obtains a simple 
relation between the critical density eo,c and the dimension 
d and argues that the critical dimension is dcrit = 10 (-1-1 
if wc include time). Our analytical approach (for a cosmo- 
logical constant A = 0) is more precise and shows that the 
critical dimension corresponding to the self-gravitating ra- 
diation has the noninteger value dc = 9.96404372... very 
close to 10 (a minor mistake in calculation was made in 
the first version of this manuscript and led to a slightly 
different value for the critical dimension. This mistake 
was pointed out to me by V. Vaganov, who performed, 
in a new Appendix B of his paper [arXiv : 0707 . 0864^4], 
a phase plane analysis of the Oppcnheimcr-Volkoff equa- 
tion with A = 0. This provides an alternative derivation 
of the asymptotic results obtained in our Sect. 15. 4p . Our 
approach also provides a thorough description of the struc- 
ture and stability of relativistic stars with a linear equa- 
tion of state p = qe for any q E [0, 1] and d > 2. Note 
that the existence of a critical dimension above which the 
oscillations of the mass-central density relation disappear 
had been noted previously in Sire & Chavanis (2002) for 
Newtonian isothermal spheres (corresponding to g — *■ 0). 
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In that case, dcrit = 10 exactly. The present study extends 
this work in general relativity. 

Acknowledgements. I am grateful to Kalyana Rama for 
pointing out his work after a first version of this paper was 
placed on arXiv:0707.2292,, 

Appendix A: Some elements of black hole 
thermodynamics 

In this appendix, we recall elementary notions of black 
holes thermodynamics in order to facilitate the compari- 
son with the results obtained in this paper. 

In the seventies, different works (Christodoulou 1970, 
Penrose &: Floyd 1971, Hawking 1971) have shown that the 
area A = AttB? of a black hole (more precisely the area 
of its event horizon) cannot decrease. Noting the analogy 
with the second law of thermodynamics, these results led 
Bekenstein (1973) to conjecture that black holes have an 
entropy proportional to their area A = AttR^ . On the basis 
of dimensional analysis, he obtained 



Sbh = A "5j ; 

-L n 



(A.l) 



where A is a dimensionless constant. If black holes have 
entropy and energy, they must possess a temperature. 
Hawking (1975) showed that black holes emit thermal ra- 
diation at a temperature 



SttGM 

Considering a Schwarzschild black hole for which 
Rc^ 



M 



2G 



(A.2) 



(A.3) 



and writing the first law of thermodynamics in the form 



d(Mc2) = 
one gets 



TdS 



BH, 



32ttXGM 



(A.4) 



(A.5) 



Comparing this expression with the expression of the 
Hawking temperature (jA.2[) it is found that A = 1/4. This 
leads to the following expression of the Black Hole (or 
Bckcnstein-Hawking) entropy H 



Sbh — 



kgA 



(A.6) 



On the other hand, using Eq. (jA.3|) . the black hole tem- 
perature (jA.5[) can be written 



kuT- 



ch 

4^' 



(A.7) 



Coincidentally, the initials are the same. 



emphasising the scaling T ^ l/R with the radius. On 
the other hand, the energy E = Mc^ of the black hole is 
related to its temperature by 



E = 



fiC' 



SnGkBT 

Therefore, black holes have negative specific heats 
dE hc^ 



C 



dT 



SnGkBT^ 



< 0. 



(A.; 



(A.9) 



Appendix B: The Tolman relation 

Global thermodynamic equilibrium requires that the red- 
shifted temperature e'^^''-'/^T(r) is uniform throughout the 
medium. This is called the Tolman (1934) relation 



e'^('')/2r(r) = To, 



(B.l) 



where Tq is a constant. Since i/(r) for r — > -|-oo, we 
conclude that Tq is the temperature measured by an ob- 
server at infinity. Let us check that this relation is satisfied 
by our equations. From Eq. (|14p . we find that 



T(r) cx e(r)'?/('?+i). 



(B.2) 



On the other hand, according to Eq. (102) of Paper I, the 
Einstein equations for a spherically distribution of matter 
give 



dP I, r.^ dy 
-d^--2^'^'^^Tr- 



(B.3) 



For the linear equation of state (O, this relation can be 
easily integrated into 



(B.4) 



e(r) = Ae i ^ , 



where A is a constant. Comparing Eqs. (jB.2p and (|B.4[) . 
we obtain the Tolman relation (|B.1[) . Now, applying this 
relation at the boundary r = R oi the system, we obtain 



T{R)=Tq. 



(B.5) 



On the other hand, according to Eq. (109) of Paper I, we 
have 



v{R) = ln(l - 2GM/Rc^). 



(B.6) 



Therefore, the temperature at infinity is related to the 
temperature at the surface of the star by 



To = r(i?)Wi- 



2GM 
Rc^ ■ 



(B.7) 



This relation shows that the thermodynamical tempera- 
ture defined by Eq. (|89p coincides with the temperature 
measured by an observer at infinity. 
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Appendix C: The Newtonian limit in a 
d-dimensional universe 

In this appendix, we briefly discuss the Newtonian limit of 
the Einstein equations in a ((i+l)-dimensional space-time. 
Considering first the general relativistic Emden equations 
(|147[) - (|148p . the Newtonian hmit corresponds to g ^ 0. In 
that hmit Eqs. pl7|) - pl8)) reduce to 



id -2) 



M(0 



and 



and they combine to give 

1 d f ^ct-i'^ 
di 



dM 
~dr 



(C.l) 



^''-1 d^ 

This is the familiar Emden equation (Chandrasekhar 
1942) with an additional factor d — 2. For d > 2, we can 
rescale the parameters so as to recover the Emden equa- 
tion exactly. For d = 2, we see that there is no Newtonian 
limit in c? = 2. This can also be seen, more fundamentally, 
at the level of the Oppenhcimer-Volkoff equations p35p - 
(jl36p . The Newtonian limit corresponds to c ^ +oo. In 
that limit, Eq. (|135p reduces to 

dP 87r(d-2) GMjr) 

dr ~ {d-\)s/ r'^-i ■ ^^-^^ 
This is to be compared with the classical condition of hy- 
drostatic equilibrium 

dP GNewtonM{r) 

= -P 



dr ^ ■ (^-^^ 

We find that the gravitational constants are related to 
each other by 

_ (rf-2)87r 

^Newton — / , i\ri ^Einstein- V^-^J 
[a — i)i>d 

They only coincide in d = 3. In d = 1. the Newtonian 
gravitational constant is infinite and it vanishes in d = 2. 
Therefore, there is apparently no Newtonian limit in one- 
and two-dimensional gravity. 

Appendix D: Dynamical stability analysis 

In this appendix, we consider the formal nonlinear dy- 
namical stability of a box-confined relativistic star with 
a linear equation of state ([8]). Specifically, we study the 
maximization problem 

Max {N[e\ \ M[e] = M fixed}, (D.l) 

and show that it provides the same condition of stability 
as the condition of linear dynamical stability studied in 
Paper I by considering the growth rate of a solution of the 
linearized Einstein equations. Therefore, linear and non- 
linear dynamical stability coincide. Stability is lost when 
the mass-central density profile A/(eo) presents an ex- 
tremum. The solutions on the series of equilibra A/(eo) are 
nonlinearly dynamically stable before the turning point of 
mass, and they become linearly dynamically unstable af- 
ter the turning point of mass. This is similar to the case 
of barotropic stars described by the Euler-Poisson system 
in Newtonian gravity (Chavanis 2002a, 2002c, 2006). 



D.l. The first-order variations 

For a linear equation of state, using Eqs. ([3]) and (fTO|) . the 
baryon number can be expressed in terms of the energy 
density in the form 



e 9+1 



1 - 



2GM{r) 



-1/2 



Airr'^dr. (D.2) 



Therefore, the first-order variations of baryon number and 
mass are 



(C.2) 5N^{^ 



1/7 



K, 



R 



'0 

1 Se 



1 - 



2GM{r) 



1-1/2 



GSMjr) 



q+l € 1 _ 2GAf(r) 



47rr dr. 



47r r" 
SM = — Ser^dr. 
c Jo 



(D.3) 



(D.4) 



Substituting SM{r) = ^ Ser'^dr' in Eq. ((D3|) . inter- 
changing the order of the integrals, and writing the first- 
order condition as Eq. ([5]), we obtain 



R 



1 



Airr'^drSel e 9+1 
\ q+l 

R 

1 

— 1 



2GM{r) 



2GM{r') 



l„2 



r'c 



-3/2 



-1/2 



Anr'dr' 



(D.5) 



This condition must be true for all variations, implying 



kV^'' 1 



-e 9+1 



2GM{r) 



_^ 2GM(r') 



-1 -3/2 



-1/2 



Anr'dr' 



(D.6) 



If we take r = i?, wc find that the value of fi is given by 



qj c2 g+1 



e{Ry 



1 - 



2GM 



Rc^ 



-1/2 



(D.7) 



On the other hand, if we take the derivative of Eq. ()D.6P 
with respect to r, we obtain, after simplification, 



2GM{r) \ de I q+l \ GM{r) AnGq 



dr 



c? q 



(D.i 



which is the Oppenheimer-Volkoff equation ^ with a lin- 
ear equation of state P = qe. This is a particular case 
of the general result given by Weinberg (1972) which is 
valid for an arbitrary equation of state provided that the 
perturbations are adiabatic. 
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D.2. The second-order variations 

We now turn to the more complicated second-order vari- 
ations. Writing 



In conclusion 
G 



2iq+l)Jo 



2qr 1 &\^2^^, 
q + 1 e dr J 



(D.9) 



dr 

the second variations of the baryon number are given by 



871- (g- 



2g + l 
9+1 



G 



R 



9+1 
+ 2^ 



dr- 



1 



1 - 



2GM{r) 



1 -3/2 



dr^ 



dr 



1 fie 



AttG 



(g + l)er 



e ,+1/ 



dr 



(D.15) 



drAi: 



2GM{r) 



-I -5/2 



47rr2 2(g -h 1)2 



2GM{r) 



rc^ 



rc^ 

^1/2 



To determine the sign of , we are led to consider the 



j2 eigenvalue equation 



29+1 / df 

e — { -f- 

dr 



in (q + 1)2 r' 



(D.IO): 



where we have omitted, for brevity, the positive term 
{q/Ky^'^ in factor on the r.h.s. This is the sum of three 
integrals that will be denoted /i, I2, and I^. Integrating 
by parts the first integral, we get 



dr^ dr 
G 1 



1 de 
e dr 



AttG 



(q + l)er 



2iq+l) r2 



3A/2 



e 9+1 1 



2qr 1 de \ 
q + 1 e dr J 



/ = A/, (D.16) 



-5/2 




with the boundary conditions /(O) = f{a) = (see Paper 
I). Introducing the variables defined in Sect. 12. 3( the fore- 
going equation can be rewritten as 



2{q + l)dr 



Expanding the derivative, we obtain 



£9 + 1 



fdr. 



de 



di) 



di 



1 + 



(D.ll) 



2q£, di! 

q + Tli^ 



i^.-^.2„-A/2 



2^G2'° 



A/, 



(D.17) 



h + h 



G 



4(9+1)70 r 



1 



2GM(r) 



(2AttG 



2GM{r) 



1 



rc^ 
2GM{r 



rc 



,2 



rc^ 

1 

qer 

-1 



-3/2 

_ g 

e 9+1 
2qr 1 de 



g + 1 e dr 
- 1 }.fdr. 



with /(O) = /(a) = 0. This equation determines a discrete 
set of eigenvalues Ai > A2 > A3, ... If all the eigenvalues 
A are negative, then < and the configuration is a 
maximum of N at fixed mass. If at least one eigenvalue is 
positive, the configuration is an unstable saddle point. We 
therefore need to determine the point of marginal stability 
(D.12) A = in the series of equilibria. It is obtained by solving 
the differential equation 



Using the Oppenheimer-Volkoff equation (jP.Sp . we find 
after simplification 



d^/j 



G 



4(9 + 1)70 r 



1 



1 



2GM(r) 



rc 



-3/2 



e 9+1 



de 

2q^ d-)p 
9+1 dC 



df 



d^ 
/ = 0, 



Aqr 1 de 



9 



1 e dr 



-- /^dr. (D.13) 



(D.18) 



with /(O) = f{a) = 0. The same equation was obtained 
On the other hand, integrating by parts the third integral P^per I [see Eq. (167)] by determining the condition of 
in Eq. (Imoll and using the Oppenheimer-Volkoff equation marginal linear dynamical stability. Therefore, the thresh- 



(|D.8[) . we obtain after simplification 



h = 



Sn (g + l)2 



drf^ 



1 



2GM{r) 



xe 



2GM{r 



dr^ dr 

-1 



-1/2 

1 de 

e dr 



AttG , 



l)er 



(D.14) 



olds of linear and nonlinear dynamical stability coincide. 
Equation (jP.lSp was solved in Paper I and led to the re- 
lation (1-171) identical to (1-148). Therefore, the points 
where A = correspond to extrema of mass in the series 
of equilibria M(a). A first eigenvalue Ai becomes posi- 
tive at the first mass peak (implying instability) , and new 
modes of stability are lost at subsequent extrema. At the 
i-th. extremum, we have Ai > A2 > ... > A^ = > A^+i > 
Ai+2 > .... 
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